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ABSTRACT 


This  study  reviews  the  history  and  development  of  marine 
electric  propulsion  drives,  the  types  of  electric  propul- 
sion, and  the  inherent  losses  which  occur  within  the  synch- 
ronous AC  machines  typically  used  for  high-power  propulsion 
systems. 

A  thorough  review  of  the  literature  pertaining  to  heat 
transfer    in      electrical   machinery   is    made.  In    particular, 

the  use  of  liguid  cooling  in  various  flow  configurations, 
including  buoyancy- driven  t hermosyphcns  and  two-phase  ther- 
losyphcns,    is  analyzed. 

Forced- liquid  cooling  is  feasible,  but  the  reguired 
rotating  seals  are  a  problem  in  reliability.  Closed-loop 
thermosyphcn  cooling  appears  feasible  at  high  rotational 
speeds,  athough  a  secondary  heat  exchange  through  the  shaft 
is  reguired.  Closed,  two-phase  thermosy phons  and  heat  pipes 
are  also  feasible,  but  reguire  forced>-air  circulation  for 
heat  rejection  tc  the  ambient.  Since  all  of  these  concepts 
deserve  additional  attention,  areas  for  further  research  and 
development    are    recommended . 
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*•  IHXBQ22CTION 

The  marine  engineer  and  naval  architect  are  charged  with 
creating  a  design  of  a  vessel  from  a  broad  set  of  require- 
ments for  the  ship  system  prepared  by  the  customer.  This 
requires  a  very  structured  procedure  for  preparing  th0 
design-based  comparisions  and  economic  trade-off  studies  at 
each  and  every  stage  cf  the  development  of  the  design.  The 
end  result  is  a  complete  set  of  design  requirements  from 
which  to  construct  the  final  product  (this  result  should 
meet    the   original  design   requirements!)  . 

This  structured  procedure  has  the  form  of  the  well-known 
"design  spiral"  [Ref.  1],  which  is  shown  in  Figure  1.1  for 
the  propulsion  machinery.  This  procedure  encompasses  the 
total    ship    design,    including    (but   net    limited   to) : 

"£.=11   Propulsion    System 

Shaft    horsepower 

Propeller    RPM 

Specific  fuel  consumption   and   bunker    capacity 

Space    and    weight   objectives 

Adaptability  to   ship   configuration 
Auxiliary   Ship    Systems 

Power    and    lighting 

Steam-galley,   deck,    and    heating   systems 

Heating,    ventilating,    and   air   conditioning 

Firefighting  and   ballasting 

Fresh    water 
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Hull   Et} gin .  e  e ting    Systems 

&  n  c  ho  r    ha  n  d  I  ir.  g 

Steering  engine  and   bridge   telemetering    control 

Cargo    handling  gear 

Crane    systems 
Electronic    a  r.d    Navigation    Systems 

Communication,    exterior   and   interior 

Radar 

Lcran,    Omega,   etc.  ,    navigation    aids 

Military  electronics,  sensors,  command  and  control 
systems,  weapons  directors,  tactical  data  systems, 
and  electronic   countermeasures. 

In  order  to  design  the  vessel  within  economic  constraints, 
or  to  use  the  life  cycle  cost  as  the  function  to  optimize, 
it  is  imperative  to  define  initially  the  devices  used  tc 
provide  the  power  for  the  vessel.  They  must  provide  most  for 
the  least,  in  terms  of  the  fuel  that  they  expend  and  the 
volume    and    weight    that   they    possess. 

The  selection  of  the  main  propulsion  plant  requires 
matching  the  power  of  the  generating  device  with  a  transmis- 
sion (in  most  cases)  ,  the  propulsor,  other  ship  systems  and 
the  hull.  If  one  limits  the  choices  of  propulsors  tc  fixed- 
pitch  propellers  and  controllable  and  reversible-pitch 
propellers,  which  are  currently  realistic  for  large  vessels, 
the  number  of  possible  permutations  remains  quite  sizable, 
as  shewn  in  Figure  1.2.  The  horsepower  and  RPM  requirements 
of  modern  warships  require  the  use  of  a  transmission  to 
couple  the  high  RPM  of  an  economical  prime  mover  with  the 
low      RPM      required      at   the      propulsor      for      high      propulsive 
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Figure   1.1        Preliminary   pesign   Spiral 
for    Propulsion   Machinery. 


efficiency.        This    transmission   may      be    either    mechanical   or 
electrical.  The        mechanical      reduction      gear        has      been 

exploited  in  the  past  to  minimize  volume/space  requirements. 
when  coupled  to  a  reversible  prime  mov=r,  such  as  a  steam 
turbine,  it  has  teen  proven  to  be  quite  acceptable. 
However,  the  electrical  transmission  has  several  attractive 
features,  such  as  ease  of  propulsor  speed  and  direction 
control,  flexibility  in  the  the  number  and  location  of 
controls,  freedom  cf  installation  and  machinery  layout, 
ccuplir.a  of  various  prime  movers  independently  to  propulsor 
drives  (as  well  as  ether  auxiliary  applications)  ,  and  the 
ease  of  auxiliary  thruster  installation  for  maneuvering 
control. 

The    choice    of   the   propulsion^    plant    includes    many    impor- 
tant   factors,    such   as: 

1 .  Reliability 

2.  Maintainability 

3.  Space    and    arrangement   requirements 

4.  Weight    requirements 

5.  Type   of   fuel  required 
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6.  Fuel   consumption 

7.  Fractional    power  and  transient   performance 

8.  Interrelations    with   auxiliaries 

9.  Reversing    capability 

10.  Operating    personnel 

11.  Faiing    limitations 

12.  Costs 
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Figure    1.2        Propulsion   Plant   Alternatives. 

These  factors  have  gone  through  many  cyclic  changes 
based  on  the  fluctuating  economy.  The  recent  emphasis  on 
fuel  costs  and  availability  has  driven  the  fuel 
con  sumption /prime  mover  efficiency  to  a  much  higher  level  of 
concern  than  it  was  previously  given.  The  cost  of  acguisi- 
tion  of  ne v  vessels  has  skyrocketed,  as  well  as  the  costs  of 
maintenance  and  repair.  These  escalating  costs  are  a  chal- 
lenge to  the  engineer  of  today.  The  economic  considerations 
of  the  initial  cost  of  a  plant,  the  cost  of  maintenance  and 
repair,        and      the        cost      to      operate      a        plant      from      its 


17 


construction  to  its  re-irement,  also  coined  "life  cycle" 
cost,  is  the  driving  force  that  controls  the  design  of  the 
■odern   warship. 

Of  the  list  of  prime  movers  shown  in  Figure  1.2,  ths 
life-cycle  cost  of  the  marine  gas  turbine  is  one  of  the 
lowest.  The  economy  of  operation  is  maximized  within  a 
rather  narrow  range  cf  RPM's.  This  has  been  tolerated  in 
the  current  employment  of  this  device,  since  the  cost  of 
utilizing  an  electric  transmission  to  accomplish  the  goal  of 
operating  the  prime  mover  within  its  maximum  economy  range 
did  not  offset  the  less  of  efficiency  experienced  by  the  gas 
turbine.  This  was  due  to  the  technology  that  existed  within 
the  field  cf  electrical  controls  at  the  xime.  The  advent  of 
electronic  frequency  control  of  the  output  of  the  generator, 
independent  of  RPM  cr  poles,  and  the  frequency  control  of 
the  synchronous  motor,  independent  of  its  supply  frequency, 
makes  the  use  of  the  electrical  transmission  much  more 
attractive.  Although  the  initial  cost  will  be  higher,  the" 
cos*  of  operation,  maintenance,  and  repair  of  the  system 
will  rapidly  offset  the  initial  cost  differential  [Ref.  2], 
[Ref.    3],      and    [Ref.    4].  An   analysis   of    the      economic   and 

engineering  considerations  involved  with  this  mechanical  vs. 
electrical  drive  selection  was  made  by  Eric  Gott  and  S.  0. 
Svensson   [Ref.    5]      in    1974;  they   clearly      illustrated   the 

nature  cf  the  problem  and  possible  solutions.  In  ^heir 
paper,  it  was  evident  that  the  selection  of  some  of  the 
possibly  mere  appropriate  forms  of  electric  drives  was 
discouraged  due  to  the  economics  of  the  choice.  High-power 
frequency-controlling  equipment  was  very  expensive  and 
bulky.  Recent  advances  in  high-power,  solid-state  equipment 
have  now  again  promoted  these  alternatives  into  viable 
contention. 
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A.       HISTORICAL    DEVELOPMENT 

From  the  aarly  attempts  by  the  Russians  in  the  1890 »s  to 
power  small  passenger  launches,  and  the  British  ar.-1.  French 
experiments  with  small  submersibles ,  to  the  present,  the  use 
of  elec-ric  propulsion  has  made  a  slow  evolution.  The  U.S. 
Navy  began  with  the  Jupiter  in  1911.  It  was  a  5500  hp/shaft 
vessel  equipped  with  wound  rotor  induction  motors,  which 
fulfilled  a  30-year  lifespan  (terminated  by  surface  warfare 
activity  in  194  3,  while  serving  as  an  aircraft  carrier 
(Largely)).  This  led  to  several  major  programs,  summarized 
belc  w  : 

*  AC  drive  installed  on  five  battle  ships  of 
30,000+hp  (1918-1922) 

*  AC  drive  applied  to  carriers  Saratoga  and  Lexington  of 
200,000hp    (1920) 

*  Synchronous  ac  drive  introduced  on  four  Coast  Guard 
cutters    (19  19) 

*  Synchronous  ac  drive  in  the  liner  Normandie  at 
160,000hp    (1935) 

*  Synchronous  ac  drive  and  full  dynamic  braking  in 
cutters     (19  40) 

In  the  U.S.  Navy  Bureau  of  Ships  Manual  (BUSHIPS 
250-660-2),  dated  1  March  1945,  the  text  states  that  the 
Navy's  inventory  of  ships  that  used  electric  drive  (in  addi- 
tion to  above)  was  "....about  two  hundred  destroyer  escorts 
and  about  one  hundred  other  ships  including  tankers,  cargo 
transports,  troop  transports,  and  store  ships."  All  of  these 
applications  were  with  steam  turbines  (which  were  speed- 
adjustable  and  reversible)  and  eventually  this  design  lest 
favor  to  the  less  expensive  mechanical  reduction  gear 
systems    which      were    somewhat    lighter.         The    weight      of    these 
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original  electric  drive  systems  was  due  to  the 
frequency  relationship  of  rem  to  the  number  of  poles.  This 
led  *:c  the  basic  configuration  of  two  pole  generators 
(3600RPH)  and        twenty-sight      to        eighty        pole        motors 

(256-9CREM)  .  The  large  number  of  poles  required  in  motor  the 
for  spe  =  d  redaction  caused  the  motors  to  be  extremely  large 
and  heavy.  The  advent  of  the  use  of  synchronous  motors 
reduced  the  weight  slightly  since  the  motor  experienced  much 
fewer  slip  losses.  A  good  historical  background  on  the 
early  stages  of  electric  propulsion  is  found  in  the  book 
written    by   Commander    S.H.    Robinson,    USN   [Ref.    6]. 

Current  applications  in    the   U.S.    Navy    are   shown    in 
Table    I. 


TABLE   I 
U.S.    Navy   Electric   Drive    Ships 

SHP 

15,000 

11 ,520 

10,000 

1  ,600 


Type 

# 

Ships 

AS 

2 

*S 

2 

TARC 

1 

TAGOS 

1 

# 

Motors 

1 

8 

2* 

2** 

*      also    has    4-12  00    hp  thrusters 
**   also   has    1-550   hp   thruster 


E.       ELECTRIC   PROPOLSICB    CONCEPTS 

An  area  of  current  research  in  the  field  of  rotating 
electric  equipment  that  shows  considerable  promise  with 
regard  to  electric  propulsion  is  the  use  of  acyclic  (homo- 
polar)  superconductive  DC  devices.  These  devices  produce 
exceptionally  high  torque  and  horsepower  in  a  very  small 
volume      and   weight.           This    is      of   great      value    in      advanced 
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hullforn   propulsion.  A   large   amount      of    research      in   this 

area  to  determine  configuration,  peak  casualty  potential, 
maintainability  of  the  required  cryogenic  equipment,  as  well 
as  the  adaptibility  of  the  systems  for  this  cryogenic  appa- 
ratus must  be  performed  before  actual  implementation  of  this 
configuration   can   occur      in    a    specific   design.  Thus,      for 

rapid  implementation  in  a  short-term  project  (3-5yrs) ,  the 
choices  remaining  in  electric  propulsion  to  be  considered 
are  direct-current  (DC)  and  alternating-current  (AC) 
devices.  The  category  of  AC  devices  is  further  subdivided 
into  inductive  and  synchronous  devices.  All  three  devices 
have  teen  used  for  trarine  propulsion  and  each  has  its  own 
inherent   advantages    and   disadvantages. 

1 •      £C    Machinery 

The  DC  electric  drive  system  has  been  used  exten- 
sively in  the  past,  primarily  in  low  to  medium  power  ranges 
(1000  -  6000  hp/shaft)  with  extremes  of  400  hp/shaft  and 
19,600      hp/shaft      noted.  This      amployment      of      DC      drives 

normally  •  utilizes    1-4     prime    movers      powering     one   or      more 
motors    per      shaft    (duplicating   the      system    if    more      than   one 
shaft   e  xists)  . 
Primary    advantages    of  this    system   include: 

1  •  IH§i§  2f  Control.  Control  being  effected  by  varying 
the  generator  voltage  through  field  control-a  process 
that   lends    itself   well    to   remote   location (s) . 

2«  Kaltiple  Control  Stations.  The  ease  of  which  the 
generator  voltage  and  polarity  of  the  motor  connections 
(either  field  cr_exciter)  are  varied  allows  multiple 
control  stations  to  be  designed  for  flexible  ship 
control  in    tight   maneuvering   situations. 

3.  adaptability  to  varinjj  propeller- hull  characteris- 
tics. The         propeller-hull      characteristics        change 
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dramatically  when  the  vessel  is  towing  or  ice-breavir.q. 
The  ability  to  match  the  prime  aov=>r  (running  at 
optimum  speed)  to  the  oropeiler  throuah  the  prooulsior 
motor  and  still  provide  maximum  torque  is  a  specific 
advantage  of  this  system,  without  excess  engine  or 
electrical  capacity,  through  a  wide  range  of  propeller 
FP3s. 
Primary    disadvantages   include: 

1 •  limited  Power  Range .  The  top-end  range  of  approxi- 
mately 10,000  hp/shaft  limits  the  vessel  size  -hat  can 
utilize  this  system. 

2-  Generator  Speed  Limitations-  The  maximum  speed 
that  the  generator  can  be  operated,  due  to  problems 
primarily  associated  with  commutation,  is  typically 
limited  to  under  1000  BPM.  This  limits  the  selection 
of  the  prime  mover  if  one  does  not  utilize  mechanical 
reduction  gearing  to  match  the  characteristics.  Figure 
1.3   shews    the   normal  relationship  of   RPM   vs.    KW. 

3-  Weight  and  Size.  The  weight  of  the  DC  system 
increases  dramatically  with  power  rating  *s  shown  in 
Figure  1.4  (this  figure  shows  weight  in  pounds/SHP)  ; 
note  that  total  weight  becomes  almost  linear  with  SHP 
at  higher  SHP  values.  This  precludes  high  power  appli- 
cations of  conventional  DC  systems  above  10,000 
hp/shaft  (acyclic-homo  polar  devices  are  excluded  from 
this  analysis) .  The  physical  size  of  the  device  also 
has  this  dramatic  functional  relationship  to  rating  as 
shown    in   Figure    1.5. 

From  the  above  discussion  of  advantages  and  disad- 
vantages, the  consideration  of  DC  electric  propulsion  for 
use  in  a  warship  design  that  would  require  on  the  order  of 
40,000   hp    and    utilize   gas   turbines   in    their    most    efficient 
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Figure    1. 4 


Variation   of   DC   Generator   Weight 
with   Bating. 


operating  range    is    net  realistic.        Therefore,   the   remainder 
cf      this      papsr    will      only      consider      the    AC      systems.         The 
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Figure  1.5    Variation  of  DC  Generator  Size 
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inherent  losses  associated  with  these   DC  systems  though  are 

l 

essentially  the  same  as  those  in  the  AC  systems.  Figures 
1.1  through  1.5  are  from  "Marine  Engineering",  R.  L. 
Harrington,  editor. 

2 .   AC  machinery 

The  devices  included  in  the  conventional  AC 
machinery  area  fall  into  two  categories:  induction  and 
synchronous.  The  principal  difference  between  the  two 
sys-.ems  is  in  the  motor  application.  The  induction  motor 
system  provides  exceptional  torque  at  low  RPM  and  maximum 
slip  (200%).  This  high  torque  response  at  low  rpm  makes  the 
induction  motor  system  preferable  in  some  applications.  For 
r.aval  applications,  this  system  is  not  as  efficient  as  the 
synchronous  system  by  a  few  critical  percentage  points.   The 
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induction  motor  has  losses  that  are  a  function  of  the  slip 
incurred  between  the  rotor  and  the  stator  field.  This  slip 
is  related  to  the  RPH  and  lead.  The  remaining  losses  in  the 
inducticr.  motor  are  the  same  as  the  synchronous  motor  in 
general.  Another      critical      difference      between      the      two 

systems  is  the  power  factor  difference.  In  the  induction 
motor,  the  rotor  current  is  induced  from  the  stater  field. 
This  induced  field  *:hax  is  created  as  a  result  of  this 
current  is  the  mechanism  from  which  the  mechanical  rotation 
is  derived.  Since  this  current  is  induced  from  the  stator 
(supplied  current)  ,  it  causes  the  current  load  seen  by  the 
supplying  dsvice  to  "lag"  the  supplying  voltage  and  there- 
fore cause  the  power  factor  (a  measure  of  the  amount  of  lag 
or  lead)  4-o  be  less  than  unity.  A  power  factor  of  unity 
indicates  that  the  load  current  and  the  load  voltage  are  in 
phase  and  maximum  real  energy  transfer  is  occurring.  The 
induction  motor  is  normally  operated  with  a  lagging  power 
factor  of  approximately  80/?  [Ref.  6]  and  an  overall  effi- 
ciency of  approximately  90%.  In  the  synchronous  system,  the 
power  factor  is  nornally  unity  and  is  controllable  .  by  the 
external' motor  rotor  excitation  which  results  in  the  synch- 
ronous system's  power  factor  being  adjustable  and  a  unity 
power  factor  being  easily  obtained.  Note  here  also  that  the 
power  factor  in  this  systam  can  be  made  to  actually  "lead" 
and  thereby  correct  a  system  "lag"  to  other  loads  that  might 
exist  in  the  fully  integrated  system.  The  advantages  of  AC 
systems    in    general    are: 

1.  High  Efficiency.  Losses  in  the  synchronous  system 
ara  approximately  6%  (including  excitation  losses)  as 
compared  to  the  induction  motor  system  losses  tha-1: 
range    from    approximately    8   to    10^. 

2.  Flexibility  of  Installation.  The  direct  drive 
motor   can    ba   installed    in      the    vessel    to    minimize    shaft 
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length  to  the  propeller  (sines  physically  smaller  than 
DC  motors)  and  the  generator  can  be  installed  without 
reduction,  gearing  and  operated  directly  coupled  to  the 
prime    mover. 

3«  Dual  Use  of  Propulsion  Power.  The  AC  system  allows 
for  the  power  from  propulsion  units  tc  be  used  for 
other  functions  when  net  required  for  propulsion;  this 
leads    to  the   creation    of   a   fully-integrated   system. 

u»  Available  in  large  Power  Ratings.  The  conventional 
AC  systems  can  be  built  in  any  power  rating  reguirec 
with  a  reasonable  upper  limit  of  approximately  60,000 
hp.  The  larger  system's  size  is  being  reduced  by 
raising  the  system  voltage  (which  is  not  limited  by  the 
commutation  reguirements  of  DC  systems).  The  system 
voltages  typically  in  use  are  from  2,300  to  about  7,500 
volts.  These  voltages  are  easily  transformed  for  other 
uses   in  addition  to   propulsion. 

The  selection  of  the  main  propulsion  system  there- 
fore is  not  simple.  The  design  environment  must  be  specified 
and  the  prime  mover  is  selected.  Once  these  parameters  are 
defined,  the  appropriate  selection  of  transmission,  control, 
and  prcpulscr  can  be  chosen.  For  the  purpose  of  the  propul- 
sion cf  a  vessel  of  approximately  8,000-ton  displacement  and 
the  mission  of  a  normal  warship  of  the  destroyer/cruiser 
variety  using  the  marine  gas  turbine  as  the  prime  mover  and 
life-cycle  costs  the  design  criteria,  the  selection  of  the 
AC  synchronous  drive  (possibly  fully  integrated)  appears  as 
a  very  appropriate  choice.  Some  of  the  design  configura- 
tions possible  within  this  choice,  and  the  associated 
performance  of  each,  are  shown  in  Figure  1.6  as  presented  in 
[Ref.  7].  Additional  possible  configurations  are  shown  in 
Figures  1.7  and  1.8  .  The  end  result  of  the  implementation 
c4-    these      technologies   would      result    in      substantial    savings 
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ever  the  current  production  configuration  that  utilizes  the 
gas  ^.urbine  main  propulsion  in  the  life-cycle  analysis.  Such 
a  comparison  is  shown  in  Figure  1.9,  also  from  the  NAVSEA 
report  [Hef.   7], 

Additional  propulsion  plant  arrangements  are  given 
in  [Eef.  5]  and  the  total  number  possible  is  only  limited 
truly  by  the  designer's  imagination  and  creativity  within 
the    confines  of    the    hullform. 

3-      i52l££/Generatcr   Losses 

The  concept  of  losses  within  the  motor  and  the 
generator  are  comparative  in  magnitude  and  virtually  iden- 
tical in  nature.  The  losses  will,  therefore,  be  discussed  as 
a  grcup  which  is  applicable  to  both  the  motor  and  the  gener- 
ator. The  approach  of  an  energy  balance  is  useful  here  to 
introduce    these    losses: 

I  Energy    into      device!         JSnergy    into    device 

I  en   the  |    +     |  on    the 

I       mechanical  side      |         |    electrical    side 


I    Increase   of 
I    stored   kinetic/ 
!    potential   mech. 
I    energy   in   device 


I    Increase   of 
I    stored   energy 
I    in    the   elec.Smag, 
I    fields   of   device 


I    * 


Heat      | 
generated! 

within 
device    I 


(1.1) 


The  net  guantity  on  the  right-hand  side  of  the  equa- 
tion is  small  and  positive;  the  majority  of  the  energy  into 
the  device  is  transmitted  through  the  device  as  an  output. 
As  an  example,  when  considering  the  motor,  the  energy  into 
the  device  on  the  mechanical  side  is  negative  in  sign  and 
approximately  equal  to  the  electric  energy  into  the  device 
which  is  positive  in  sign.  This  amount  must  be  balanced  by 
the    guantity  on    the    right  side   of   the   equation. 
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The  terms  on  the  right  side  are  all  usually  small  and  may  be 
zero.  The  first  tern  represents  the  energy  associated  with 
a  scee^  chance  and  results  ^r^  th  e  in  =>^~  j  a  of  ~he  5  -v'  •~~- 
Ir.  the  steady  state  it  is  zero.  The  second  term  would  result 
fro  a  a  change  in  current  to  the  device  and  the  electric 
field  that  sxisted  would  be  a  source  of  energy  (the  induc- 
tion effect).  The  third  term  is  the  main  point  of  interest 
here.    It    is   comprised  of  the   following: 

*  Mechanical   friction    losses 

*  Air    friction    losses     (windage) 

*  Hysteresis  and  eddy-current   losses 

*  Resistance   losses   in    the   armature,    and 

*  Exciter  losses  (in  the  case  of  synchronous  machines) . 
The  object  for  the  electrical  engineer  is  to  minimize  these 
losses  tc  the  greatest  extent  possible.  Since  it  is  impos- 
sible to  eliminate  all  these  losses,  the  problem  that 
remains  is  the  removal  of  the  resulting  heat  within  the 
device  [Ref.  8].  The  size  and  weight  of  the  device  is 
largely  dependent  on  the  ability  to  remove  this  heat,  since 
the  capacity  of  the  device  is  limited  by  the  insulation  used 
in  the  device.  The  method  of  heat  removal  is  not  limited  to 
a  single  mode,  but  is  a  combination  of  many  modes  of  heat 
transfer   including: 

*  Conduction:       stator    to   casing 

rotor  to   stator    (through   air    gap) 
rotor  to    shaft 
shaft   to   casing 

shaft  to   external   device    (coupling, 
reduction   gear,    etc.) 

*  Convection    (both   forced  and   natural) : 

casing   to   ambient 

stator    tc    internal   air    (or   gas) 

rotor  to   internal   air    (or  gas) 

shaft   surface   to    internal    air    (or   gas) 
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shaft   surface   to    external    air 
*   Radiation:        external   surfaces   to   sinks 

internal   surfaces   to    casing 

These  modes  are  not  all  very  controllable.  large 
amounts  cf  research  have  been  devoted  to  minimizing  these 
thermal  resistances.  Efforts  to  improve  convection  heat 
transfer  have  found  mat  the  restrictions  on  the  rotor- 
stafor  air  ga?  have  prevented  exploitation  cf  this  mode  as 
discussed  as  early  as  in  1926  by  G.E.  Luke  [Ref.  9]  and  as 
late  as  in  1979  by  O.N.  Kostikiv,  et  al.  [Ref.  10].  The 
improvement  of  the  conduction  heat  transfer  has  been 
improved  by  minimizing  the  thermal  resistance  cf  the  insu- 
lating compounds  and  the  materials  used  in  the  construction 
cf  the  device,  including  installing  thermally-conductive 
materials    in  the   coil   ends. 

The  improvement  of  convection  heat  transfer  is  still 
being  examined  by  using  liguids  and  gases  to  cool  the  rotor 
and  the  stater  of  these  devices.  This  procedure  involves 
piping  a  liguid  through  the  stator  to  remove  heat  by  forced 
convection,  and  utilizing  a  material  with  high  thermal 
conductivity  directly  adjacent  to  the  piping  to  conduct  it 
to  the  piping  where  it  is  removed  by  the  passing  liquid. 
This  method  is  equally  applicable  to  the  rotor,  provided  the 
problem  cf  the  rotatinq  seals  can  be  overcome  in  order  to 
channel  the  liquid  from  an  external  sink  to  the  moving  rotor 
and  back  out.  The  advantage  of  doing  this  over  standard 
cooling  schemes  is  that  the  capacity  of  a  device  may  be 
increased  due  to  the  lower  internal  temperatures  that  could 
be  maintained  (or  the  physical  size  and  weight  of  the  device 
could  be  reduced  fcr  an  existing  capacity)  as  shown  in 
Pigure  1.10  for  a  40,000  hp ,  180  RPM  synchronous  motor.  The 
use  cf  a  liquid  in  a  bouyancy  driven  closed-loop,  called  a 
ther mosyphen,    is   also   possible.      The   amount   of  heat    transfer 
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in    this    Bethcd    is   dependent    or.   the   specific   hsat    capacity    of 
the    liquid    being      used,      the   temperature    rise      allowed,      the 


*    _  i 


pipe    resistance    m       .he   loop,      and    the      required    floi 

the    heat-transfer      areas   then   required   may      be   prohibitively 

large. 

The  two-phase  cooling  method  within  a  closed  loop, 
in  the  form  of  a  heat  pipe,  has  also  been  extensively  inves- 
tigated. This  method  has  an  advantage  over  the  conventional, 
liqui d-ccoling  method  in  that  it  eliminates  the  necessity 
for  rotating  seals.  Thess  seals  can  be  weak  points  in  the 
design  from  a  reliability  stand  point,  especially  if  high 
RPMs  are  involved.  The  heat -pipe  method  typically  involves  a 
confined  fluid  acting  as  a  two-phase  medium  for  heat 
transfer,  exploiting  the  phase  change  as  a  vehicle  for 
substantial  heat  transfer  in  order  to  remove  heat  from  the 
rotor  or  the  stator  or  both.  This  method  then  transfers  the 
heat  either  to  the  ambient  air  or  to  some  other  gas,  liquid 
cr  solid  heat  exchanger.  The  configuration  most  suitable 
for  use  within  a  motor/generator  on  the  order  of  40,000  hp 
would  be  either  to  ambient  air  or  to  a  fresh  water  heat 
exchanger  external  to  the  motor  in  the  case  of  stator 
cooling,  or  to  internal  forced  air  through  extended  surfaces 
in   the   case    of    rotor   cooling. 

C.       PBOBIEM    DEFINITICH 

The  purpose  of  this  thesis  is  to  analyze  various  methods 
of  cooling  motors  and  generators  of  warship  propulsion 
plants  and  to  discuss  the  economic  and  reliability  perspec- 
tives associated  with  each  system.  The  next  two  chapters 
will  discuss  the  theory  and  application  of  these  two 
methods:  liquid  cooling  of  the  rotor,  both  forced  convection 
and    the    thermosy phon ,   and   two-phase    cooling   of  the    rotor. 
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A  test  apparatus  for  evaluation  oc  va^io'i^  r z tat  ino 
ceding  schemes  has  been  designed  and  is  currently  under 
construction  at  the  David  ff.  Taylor  Naval  Ship  Research  and 
Developnent  Center,  Annapolis,  Maryland.  It  is  sketched  in 
Figure  1.11  and  will  be  used  for  the  analysis  within  this 
study.  The  results/recommendations  herein  may  thereby  be 
experimentally    evaluated. 

For  the  purpose  of  this  arbitrary  configuration,  the 
model    for  analysis    is  as   follows: 

1)  The    device  will   be  a    water-cooled    frame, 
synchronous    AC  generator,    25,118    kVA    3    3,600    RPM. 

2)  A   specific  conduction   bar    shall    be    analyzed    for    a 
typical    load    for    all    cooling    configurations.    It    has   a 
length   of    0.9144m   and   is    0.0116m    sguarewith 

a   4.763mm    diameter   hole   bored    through    the    length 
(herein    referred  to    as   the   tube   length   and    tube 
diameter) .      The   losses    per    bar    in    this    configuration 
are    approximately   50W. 
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Single-phase  cooling  has  been  extensively  investigated 
with  regard  to  rotating  references.  Work  was  stimulated 
early  by  the  printing  industry,  with  the  need  to  cool 
printing  roll  mills.  The  work  was  further  encouraged  by  the 
gas  turbine  industry  with  the  need  to  cool  high-temperature 
turbine  blades,  which  rotate  at  high  RPM,  and  by  the  elec- 
trical machine  industry,  where  the  need  to  cool  large, 
powerful  generators  became  a  necessity.  Research  in  the 
study  of  flow  and  heat  transfer  in  a  rotating  reference 
system  is  applicable  to  forced-convection  flow,  buoyant  flow 
in  thermcsyphons ,  and  the  single-phase  flows  of  closed  two- 
phase    devices. 

A  theoretical  analysis  of  flow  in  a  heated,  horizontal 
tube  without  rotation  by  Morton  [Ref.  11]  was  completed  in 
1957.  Although  his  analysis  was  limited  to  non-rotating 
flows,  it  gave  considerable  insight  into  buoyancy-induced 
secondary      flows.  His      analysis    was      limited      to      laminar 

convection  in  horizontal,  uniformly-heated  tubes  at  low 
Payleigh  numbers  (bas<=>d  on  the  temperature  gradient  along 
the  pipe  wall) .  The  work  was  an  extension  of  the  study  by 
Nusselt  who  ignored  the  gravitational  (acceleraticnal) 
effects  and,  therefore,  was  independent  of  orientation. 
This  study  was  a  numerical  solution  to  the  governing  egua- 
tions  assuming  small  rates  of  heating  and  a  notable  varia- 
tion of  density  along  the  pipe  length  giving  rise  to 
secondary    flows.  A  similar      study    was      completed    by      Kays 

£Ref.    12]      in    19  55.  Morton      studied      this    secondary      flow 

driven  by  the  buoyancy  effects  and  the  resultant  temperature 
distribution  by  successive  approximation  in  a  power  series 
(this   is      similar   to      the   flow      in   a      pipe   rotating      about   a 
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perpendicular   axis).  Seccrd-order   terms    were      required   to 

account    for      the   secondary     flow.  The    basic      equations   and 

3e  v@  1  c  c  m  e  n  t  a T*  -  or  e  s  ^  v^~  **  d  ir  RDDenfliY  ^  arifl  pre  —  he  v?ci  <? 
for  many  la^er  studies.  This  development  is  adapted  to  the 
rotating  case  by  the  substitution  of  centripetal  accelera- 
tion   for   the  gravitational    acceleration    term. 

The  configuration  of  the  cooling  loop  has  three  distinct 
regions  which  must  be  considered:  the  radial  sections  of 
tubing  (oriented  along  the  line  of  the  centripetal  accelera- 
tion) ,  the  entrance  regions  of  the  axial  sections,  and  the 
axial  sections  through  the  conductor  bars  (oriented  parallel 
to  the  axis  of  rotation).  Additionally,  it  mast  be  deter- 
mined whether  the  flow  is  laminar  or  turbulent.  The  devel- 
opment will  be  reviewed  chronologically  within  each  topic 
with  sample  calculations  made  for  pertinent  correlations  in 
Appendices    E  and  C    as  noted. 

A.       HISTOBICAL    DEVELOPMENT    OF    ONCE-THROUGH    COOLING 

In  the  196  8  publication,  "Recent  Advances  in  Heat 
Transfer"  [Ref.  13],  Kreith  made  a  survey  of  virtually  all 
the  work  to  that  date  that  had  transpired  in  the  field  of 
convection  heat  transfer  in  rotating  systems.  His  section 
en  concentric  cylinders  and  rotating  tubes  summarized  the 
works  mentioned  herein,  although  principally  evaluating  the 
reports  using  gases  as  the  fluid  media.  He  presents  the 
concept  of  utilizing  a  thermosyphon  for  the  cooling  of 
rotors,  the  theory  cf  which  is  intimately  related  to  the 
previous  work  with  single-phase  heat  transfer,  which  will  be 
discussed  later.  Kreith  also  describes  the  three  different 
sections  (radial,  entrance,  axial)  of  the  device  which  this 
study  is  addressing.  These  have  been  analyzed  individually 
both  experraentally  and  theoretically.  The  inlet  region  has 
the    swirl      effect   due      to   the      Coriolis    force,         which    is      a 
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function  of  the  Rossfcy  number,  or  the  swirl  number  (the 
reciprocal  of  the  Rossby  number);  the  radial  and  axial  tube 
sections  have,  when  a  significant  centrifugal  accelerat ior 
term  is  present,  noticeable  hydrostatic  effects  and,  addi- 
tionally, buoyant  effects  when  a  density  differential 
exists. 

1.      Ba  dial    Sections 

Extending  his  research  to  a  vertical  tube  rotating 
about  a  perpendicular  axis,  Morris  (Ref.  14]  confirmed  that 
the  approximate  solutions  must  include  second-order  terms  by 
studying  low  heating  rates  and  using  a  method  similar  to 
Morton,  but  including  tangential  terms  as  well  as  radial 
acceleration  terms.  He  identified  the  three  components  of 
convection: 

1)  Forced    convection    due   to   the      externally-generated 
pressure   gradient, 

2)  Gravitational   buoyancy   in   the   axial   direction,    and 

3)  Rotational    buoyancy   due   to   centrifugal   and   Coriolis 
forces. 

The  relative  magnitude  of  the  terms  due  to  item  3)  is  char- 
acterized by  the  Possby  number.  This  is  a  ratio  of  th<=> 
inertia  force  to  the  Coriolis  force.  Although  the  Rossby 
number  exists  in  the  velocity  and  temperature  fields  for 
solutions  through  the  second  order,  it  has  no  effect  on  the 
solutions   for  heat      transfer   or   flow   resistance.  Since   at 

the  nominal  rotation  rates  and  eccentricity  that  exist  in 
rotating  electric  machinery,  the  centrifugal  acceleration  is 
much  much  greater  than  the  gravitational  acceleration,  the 
former   will   dominate. 

In  their  1968  report  on  heat  transfer  in  rotating, 
radial,  circular  pipes  in  the  laminar  region,  Mori  and 
Hakayana  (Ref.  15]  continued  their  research  by  evaluating 
steady,         fully-developed   flow      with    analytical      techniques. 
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Using      *he    boundary- layer      concept,        they   correlated 
results      with   ex ceri mental      work   by      Trsfethen.        They 


gave 


correlations  for  the  (Su/Nu0)  ratio  as  a  function  ..  :he 
angular  velocity,  Reynolds,  Prandtl,  and  rotational  Reynolds 
numbers  for  fluids  with  Prandtl  numbers  above  and  below 
unity  (i.e.,  gases,  liquids,  glycerol,  and  liquid  metals)  . 
Their  s-udy  confirms  the  stabilizing  effect  of  the  secondary 
flow  en  the  critical  Reynolds  number  and  notes  that  laminar 
solutions  are  still  valid  when  the  Reynolds  number  is 
considerably  large     (12,600    a   n    =500    rad/s) . 

As  a  second  report  on  the  subject  of  rotating  radial 
pipes,  Mori,  et  al.  [Ref.  16]  evaluated  secondary  flow 
effects  due  zo  the  Coriolis  force  by  assuming  a  boundary 
layer  alcng  the  wall.  With  experimental  results  and  corre- 
lations from  previous  work,  they  were  able  to  show  that  the 
effect  of  secondary  flow  is  less  in  the  turbulent  region 
than  in  the  laminar  region,  although  still  on  the  order  of  a 
10*  increase  over  the  value  without  secondary  flow.  The 
correlations  they  presented  agree  with  their  previous  report 
[Ref.  15]  and  with  experimental  results  with  air;  they  are 
presented  as  =  ratio  of  Nusselt  numbers  based  on  an  exponen- 
tial function  of  the  Reynolds  and  Prandtl  numbers  and  are 
given  for  a  wide  range  of  Prandtl  numbers  including  gases 
and    liquids. 

Ito  and  Nanbu  [Ref.  17]  experimentally  measured 
friction  factors  for  flow  in  rotating,  straight  pipes  of 
circular  cross-sect icn,  and  compared  their  results  to  those 
of  Trsfethen  and  Mori  &  Nakayama.  Their  empirical  correla- 
tions for  friction  factor  ratios  for  two  ranges  of  Kt  (where 
Kt  =  J2/Re  and  J=  Q.  d2/  v  )  are  shown  in  equations  (2.1) 
and  (2.2)  for  the  turbulent  region,  Kt  >  Ktcr.  These  are  as 
follows : 


U1 


0.942    Kto.as2f     1<    Kt    <500, 


(2.1) 


■p  7f    — 


0.942     Kto-osf         Kt>500. 


(2.2) 


Eelow    a      Kt      of    1,    the   friction   factor   practically    coincides 
with    the   stationary    friction   factor   giver,    by   Blasius: 


f      =   0.316  Re- o.25 . 

0L 


(2.3) 


For    laminar   flow,      which   occurs   whan      Kt 
fcllcwina   correlation  aDolies: 


is    b^low    Ktcr,      the 


f/f      =    0.0883     Klo«2S(1    *     11.2    Kl-o.325)r 


(2.U) 


where        Kl      =        J*F.e.  Equation         (2.4)  is      valid        fc: 

2.2*102<   Kl   <10?      and      (J/Re)<0.5.  The      friction      facto: 

without   rotation   is    given   by: 


f  0     =   64/Ee 


(2.5) 


Transition    occurs   at 


Ktcr    =    1.  C7   Ji.23*  1C3 


(2.6) 


or 


Re      =1  .07    jo  .  z  3*1  o3 

cr 


(2.7) 


for    the     range    28<   J    <2*103 .  Calculation   of     the   pressure 

drop  would  be  required  in  the  detailed  calculations  for 
design  of  a  cooling  system,  especially  of  the  closed, 
rotating-loop  system. 
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2.      Axia].   Sections 

Kuo#  et  al.  [Ref.  18]  analyzed  the  heat-transfer 
characteristics  of  water  flowing  through  a  partially-filled 
pipe  which  was  rotating  about  its  own  axis  and  which  was 
maintained  at  a  constant  heat  flux.  They  also  summarized 
the  previous  pertinent  works  that  were  available  at  the  time 
the  paper  was  written.  The  paper  illustrates  hew  the  tran- 
sition from  gravitational  to  centrifugal  force  is  accompa- 
nied by  instabilities  as  previously  discussed  by  Taylor 
[Ref.  19]  and  [Ref.  20].  The  authors  calculated  a  critical 
Taylcr  number  at  which  the  Nusselt  number  changes  its  func- 
tional dependence  on  the  Reynolds  number.  The  paper  also 
shows  how  the  heat  transfer  is  accomplished  within  a  layer 
around  the  periphery  of  the  pipe.  A  correlation  for  this 
film  thickness  is  provided,  and  experimental  results  confirm 
its  validity  and  illustrate  how  the  increasing  depth  of  the 
film  in  the  rotating  pipe,  up  to  completely  full,  does  not 
enhance  the  heat  transfer  beyond  the  critical  depth  value. 
The  range  cf  the  rotational  Reynolds  numbers  used  in  the 
experimental  work  was  between  4,000  and  20,000;  the  enhance- 
ment      (Nu/Nu  )     was    approximately    up   to    240^. 

Cliver  [Ref.  21]  analyzed  natural  convection  in  a 
horizontal  tube  with  no  rotation,  following  the  work  of 
Colburn  (Ref-  22],  This  study,  while  pointing  out  several 
problems  with  the  original  research  and  basic  criticisms 
noted  by  Hartinelli  [Ref.  23],  combines  the  work  of  many 
authors  to  obtain  a  correlation  (eguation  (2.8))  for  the 
Nusselt  number  in  terms  of  the  Graetz,  Grashof,  and  Prandtl 
numbers  and  the  L/D  ratio.  The  paper  also  points  out  that 
the  Grashof  number  must  be  evaluated  using  the  temperature 
gradient  along  the  pipe  axis  instead  of  the  gradient  from 
the  pipe  wall  to  the  fluid  mean  temperature. 
Nu       (;jw/ub)0-  i*    =    1.75  {Gzm+5. 6*10"** 
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(Gr     Pr       (1/01)0.70(0,333, 

m     ra 


(2  .Si 


The  correlation  is  supposed  to  incorporate  the  effects  of 
natural  convection  along  with  forced  convection  by  vector 
addition,  tut  he  nctes  that  the  velocity  profile  of  the 
fluid  flew  is  critical  in  the  analysis  and  that  natural 
convection  in  the  case  of  horizontal  tubes  should  be  inde- 
pendent of  I/D  beyond  the  entry  region.  For  this  concept,  a 
correlation  is  also  provided  for  the  Nusselt  number  without 
the    L/D    term: 


Nu        (yw/yb)0.  1  ♦    =     1.75  {Gz 


m 


+    O.C083(Gr   Pr   )o. 75)0.333. 

m      m 


(2.9) 


The  study  of  the  Taylor  vortices  was  continued  by 
Pattenden  [fief.  24 ]  using  a  horizontal  (axial)  tube  rotating 
about  its  own  axis.  Due  to  problems  in  the  slip-ring  appa- 
ratus used  in  the  experimental  wort,  the  errors  in  the  accu- 
rate measurement  of  the  temperature  precluded  specific 
determination  of  exponents  for  the  Reynolds  number  in  the 
correlation  for  the  heat-transfer-coefficient  equation. 
Qualitative  analysis  of  the  enhancement  due  to  tube  rotation 
was  made,  noting  that  the  axial-flow  effects  were  small 
compared  to  the  rotational  effects.  Taylor  vortices  were 
neither  confirmed,  nor  denied;  the  flow  was  sufficiently 
turbulent  to  mask  the  effects  that  the  vortices  may  have 
had.      Rotation    rat»s    were   varied    to    a    maximum   of    4,000    RPM. 

Mori,  et  al.  [Ref.  25]  investigated  further  into 
the  effects  of  buoyancy  on  forced-convective  heat  transfer 
in  horizontally-oriented,  heated  tubas  with  laminar  flow. 
His  study  extended  the  work  of  Graatz  and  Nusselt,  and  went 
beyond  the  limitation  of  small  heat  fluxes  such  as  those 
studied   by   Colburn,    Beldams,    Martinelli,    at   al.,    Jackson,    et 
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al . ,  ana  Oliver.  The  majority  of  the  sxperimental  work  by 
the  anthers  was  performed  with  a  constant  wall  temperature 
and  large  L/D  ratio.  Their  study  involved  a  iargs  tempera- 
ture difference  between  the  wall  and  the  fluid  such  as  would 
occur    with    large      heat   flux.        The   change      in    viscosity   with 

temperature   is    approximately  covered      with    the    (p      /y      )o.i* 

w      *"  b 
term    m    equation    (2.8),    but    the    effect   of   natural    convection 

has    not    been  included  specifically.  This   effect    was    shown 

[Ref.  11]  as  a  function  of  the  ReRa  product;  it  is  appli- 
cable in  the  range  of  RePa  less  than  3,000,  but  within  this 
range  the  effect  of  natural  convection  does  not  exceed 
several  percent.  Tfcs  Rayleigh  number  is  a  function  of  the 
fourth  power  of  the  diameter;  therefore,  the  effect  that 
natural  convection  makes  in  small-diameter  channels  is 
small.  This  illustrates  that  in  the  thin  layer  near  the 
wall,  the  temperature  and  velocity  change  abruptly  and  are 
guite  different  from  that  in  the  inner  part  of  the  fluid. 
When  the  ReRa  product  is  large,  the  secondary  flow  is  strong 
and  the  velocity  distribution  in  the  axial  direction  is 
utterly  different  frcm  the  Poiseuille  flow  assumed  in  the 
previous    works.      The    Nusselt   number    is    then  calculated    by: 


Nu    =    2aq/   k(t    -t    ) 
w     m 


(2.10) 


where    the      stationary  Nusselt    number    for      constant   heat^flux 
case    is: 


Nu  0  =    U8/11. 


(2.11) 


The  experiments  were  conducted  with  air  and  fluids  with 
higher  Prandtl  numbers;  they  were  only  considering  heating 
in  a  gravitational  field  (not  an  acceleration  field  due  to 
rotation).  In  this  case,  the  layer  near  the  wall  becomes 
thinner   and   the    Nusselt    number   increases   with   the    increasing 
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ReRa.  The  local  Nusselt  number  for  the  experimental  ccrra- 
laticn    (fcr   air)     is    given  by: 

Ha    =   0.61  (HeRa)o. 2  (1+1.8/ (EeEa)  o«2)  {2,12} 

This  study  was  extended  into  the  turbulent  region,  where  the 
highest  point  of  the  temperature  distribution  was  found  to 
shift  in  the  direction  of  gravity  slightly.  In  this  case, 
the  effect  of  the  secondary  flow  was  not  as  significant  as 
in  the  case  of  laminar  flow.  The  Nusselt  numbers  calculated 
within  this  regime  were  in  good  agreement  with  those  of 
Colburn    [Ref.    22],    where: 

Nu    =    0.0204    Reo-s.  (2.  13) 

Mori  and  Nakayama  used  the  theoretical  analysis  of 
their  previous  research  [Bef.  26]  and  [Ref.  27]  on  straight 
pipes  rotating  about  a  parallel  axis  [Ref.  28].  Analyzing 
the  body  forces  driving  the  secondary  flow  caused  by  density 
differences  in  the  centrifugal  field  and  the  Coriolis  force 
with  regard  to  the  flew,  they  used  fundamental  principles  to 
characterize  the  flow  and  temperature  fields.  This  was  done 
assuming  an  effective  secondary  flow  due  to  buoyancy  and 
using  the  rotational  Reynolds  number.  This  would  indicate  a 
rapid  divergence  frcm  the  results  of  Morris  [Ref.  14]  and 
shows    that    Coriolis      effects   cannot    be    ignored.  The    paper 

gave  correction  factors  as  a  function  of  Reynolds,  rota- 
tional Reynolds,  and  rotational  Rayleigh  numbers,  and 
concluded  that  the  Coriolis  effects  diminish  with  increasing 
eccentricity  (higher  g-fields)  to  a  negligible  value  when 
the    ReRa    product  is    greater    than    10,000,000. 

In  1968,  Nakayama  [Bef.  29]  further  analyzed  a  hori- 
zontal, straight  pipe  rotating  about  a  parallel  axis.  He 
assumed,         as      previously      [Bef.    15],  that      an      effective 
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secondary  flow  in  fully-d eveloped  conditions  of  flow  and 
temperature  fields  existed,  and  included  body-force  terms 
explicitly.  The  correlations  he  obtained  for  both  the  ratio 
of  friction  factors  and  the  ratio  of  Nusselt  numbers  are 
valid  cv=r  a  large  range  of  Prandtl  numbers  for  both  liquids 
and  gases,  as  well  as  for  a  wide  range  of  Reynolds  and 
Grashcf  cumbers.  The  correlation  for  the  Nusselt  number  for 
liquids   in    fully-developed    flow   is: 

Nu    =    3e°-9    ?r°-*     {0.033    ( Re/  r 2 • 5) 1/3° 

(1  +  0.01U/(Re/   r2,5)  l/6))  (2.14) 

where        F   =   Re22/13     (GrrPr0  •  6)  ~l  2A  3 .         The      r        term    is   the 
ratio   of      inertia   force     to    body      force.         Inertia      force   is 
represented     by    Rs22/i3f         and   the      remaining     terms    in      V 
represent   the      body    force.  The   body      force    is      either   the 

Coriclis  force  in  the  case  of  the  inlet  region,  curved 
regions,  or  the  radial  arms,  or  the  centrifugal  force  in  the 
case  cf  the  axial  section.  The  numerical  result  of  his 
correlations  for  the  model  presented  in  Chapter  I  is 
included  in  Appendix  B  for  both  turbulent  and  laminar 
liquids. 

Sieawarth,  et  al.  [Ref.  30]  further  evaluated  the 
effect  that  bou  yancy-induced  secondary  flow  would  have  on 
laminar  flow  in  a  heated,  horizontal  pipe  without  rotation. 
His  work  shewed  clearly  that  secondary  flows  did  exist  due 
to  density  variations  and  that  the  heat  transfer  was 
enhanced  by  them  as  a  function  of  the  Grashof  and  Prandtl 
numbers.  This  effect  is  somewhat  negated  in  the  presence  of 
extreme  acceleration  as  that  in  the  periphery  of  a  rotor, 
but  remains  nonetheless,  provided  the  turbulence  is  below 
the    critical  Revnolds   value. 
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With  no  assumptions  as  to  the  structure  of  the  flow 
and  temperature  fields  to  simplify  the  governing  equations, 
Seeds   and   Harris   presented    a     numerical   solution    to 


of  the  rotating  cylindrical  tube  [Ref.  31]  with  laminar, 
fully-developed  flow.  Data  were  compared  to  the  theoretical 
results  for  the  case  of  air,  water,  and  glycerol.  The 
correlation  obtained  is  in  terms  of  the 

Bayleigh-P.eynolde-Prandti  product   as    is   given   by: 

Nu    =    Nu0  {0.262  (RaRePr)  o-1  73}  .  (2.15) 

The  value  for  the  test  model  is  calculated  and  presented  in 
Appendix  B. 

Stephenson  again  studied  this  parallel,  rotating- 
pipe,  heat-transfer  problem  for  fully-developed  turbulen*: 
flow.  He  compared  his  experimental  results  and  correlation 
with  the  earlier  work  of  Morris  and  woods  [Ref.  32]  and  they 
compared  favorably  within  the  turbulent  entrance  region.  He 
noted  that  the  rotational  buoyancy  was  net  a  strong  factor 
in  ^he  secondary- flow  in  the  fully-developed  region. 
Because  of  this,  his  results,  when  compared  against 
Nakayama's  results  assuming  a  strong  effect  [Ref.  29],  did 
not    compare   well.      His   correlation   is    listed   below: 

Nu    =    0.0071    Reo-88    jo. 023  (2.16) 

and  is  compared  with  previous  correlations  for  the  test 
model    in    Appendix  B. 

3 •      Entrance    Regions 

In  1969,  the  Institution  of  Mechanical  Engineers, 
Thermodynamics  and  Fluid  Mechanics  Group,  sponsored  a  sympo- 
sium, en  the  subject  of  heat  transfer  and  fluid  flow  in 
electrical    machines.  Included   in   the   presentations      was  a 
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paper   by     Davidson   [  Bef .    33],      which   described 
gene  rater-rotor    test      rig   for      evaluating   the      heat-transfer 
characteristics    of      hydrogen  cooling.  His  testing      shoved 

that  enhancement  in  the  turbulent  flow  regimes  was  net  as 
high    as   theoretically  predicted.  It   was    hypothesized  that 

♦he  reason  was  the  inability  cf  the  hydraulic-diameter 
concept  to  accurately  account  for  the  flow  conditions  in  the 
non-circular  ducts.  The    tests      were   able      to   give      a   good 

turbulent  correlation  with  theoretical  prediction  when  the 
cooling  scheme  was  modified  by  shortening  the  axial  path. 
This  effect,  although  not  explained,  could  have  been  due  to 
the  Coriolis  swirl  effect  in  the  inlet  region  enhancing  the 
heat  transfer  above  that  to  be  obtained  in  the  fully- 
develcped    region. 

The  effects  cf  the  entry  length,  especially  with  low 
Prandtl  numbers  such  as  those  of  gases,  are  very  difficult 
to  eliminate  and  the  Corio lis- effects  on  the  heat  transfer 
may      be      noted      even      at   high      Reynolds      numbers.  As      the 

Eeynolds-Rayleigh  product  increases  above  about  1,000,000, 
there  is  a  tendency  for  the  amount  of  enhancement  to 
diminish.  This  is  attributed  to  the  turbulent  effects  over- 
iding   the  secondary    flow   enhancement.  The    effect    of    rota- 

tion on  heat  transfer  is  to  enhance  it  in  the  entrance 
regions  ty  the  secondary  flow  due  to  the  Coriolis  effects 
and  in  the  fully-developed,  horizontal-tube  regions  by  buoy- 
ancy effects.  A  sitrilar  survey  of  the  technology  was  mada 
by  Petukhov  and  Polyakov  in  the  USSR  with  much  the  same 
results      [Ref.    34],  The      heat-transfer        problem      in      the 

entrance  region  of  a  tube,  where  the  Coriolis  terms  domi- 
nate, is  addressed  by  Morris  and  Woods  [Ref.  32]. 
Correlations  are  presented  for  air  and  are  applicable  for 
ether    gases     as    well.  These    correlations,         for    both      the 

laminar  and  the  turbulent  cases,  are  presented  as  functions 
of        the         product         of      the        Reynolds        number         and        the 


U9 


rotational— Reynolds      number.        The     Morris      and   Woods 


notes  that  further  work  is  required  in  the  case  of  liquids, 
although  a  similar  approach  is  valid.  Mori  ar.d  Nakavama 
also  presented  a  survey  of  the  state-of-the-art  technology 
regarding  the  heat-transfer  characteristics  of  rotating 
pipes  and  ducts  [  Ref .  35].  The  questions  of  the  entrance 
effect  and  its  length  were  emphasized  as  being  generally 
unanswered,  but  they  gave  a  general  guideline  for  deter- 
mining the  extent  of  the  entrance  region  as  approximately  20 
times  the  diameter  of  the  pipe.  They  also  verified  correla- 
tions  cf   previous   works   for    helium   and   water. 

4.      Combinations    cf     Radial,      Horizontal.        and    Entrance 
Sections 

Also  presented  in  the  1969  Institution  of  Mechanical 
Engineers  Symposium  was  a  study  by  Lambrecht  (Ref.  36], 
which  discussed  the  problems  associated  with  water  cooling 
cf  rotors.  His  paper  summarized  the  theoretical  considera- 
tions, as  well  as  reviewed  his  previous  work  and  the  work  of 
fellow  German  researchers  (Neidhoeffer  and  Ingenieure) .  He 
noted  the  superior  ccoling  properties  of  water  and  the  size 
and  weight  reduction,  along  with  the  improved  efficiency  of 
water-cccled  machines.  He   listed      values   for      the   optimum 

cooling  duct  size  based  en  the  heat-transfer  and 
electrical-loss      characteristics      for    various     gases.  The 

method  used  for  optimum  duct  calculation  was  based  on  the 
principle  that  the  pressure  drop  in  the  duct  decreases  with 
increasing  duct  size  and  electrical  losses  increase  with 
decreasing  conductor  area  (increasing  duct  size).  A  suit- 
able calculation  of  this  type  is  necessary  in  any  final 
design  for  cooling  of  these  devices  by  any  method  devised 
for    either    the    single-   or   the    two-phase   approach. 

In  1970,  Sakamoto  and  Fukui  measured  heat-  and 
iass-transf er   coefficients    for   air   and   oil,    specifically   for 


50 


use  in  ceding  rotating  electric  machines  [Bef.  37].  Noting 
that    the  roost    effective   and      economic    method    to   increase    the 

C  C  W  3  X         0  C^  P  VL^l     C  ~  3  2.  £  C^  IT  2.C  &  1         "tSLC^-^Tl^^         <-r  ^s  c;     *■  q         t   ri  ?*>  ^~  f  ?r  c:         <*">■*- 

cooling  cf  the  insulating  material  (especially  in  the  rotor 
conductors  and  drums) ,  a  gscmetry  similar  to  a  typical  rotor 
was       fabricated    and      suitable      coolants      were   tested.  The 

convective  heat  transfer  was  characterized  as  a  function  of 
the  L/d,  Re ,  3rr,  Pr  variables,  with  the  H/d  parameter  held 
constant.  When  the  shape  factor  of  eccentricity  H/d  is 
sufficiently  large,  the  centrifugal  acceleration  may  be 
assumed  to  be  a  uniform  acceleration  (  2H)  acting  on  the 
tube  axis,  which  is  similar  to  the  force  gravity  exerts  on  a 
horizontal  tube.  with  a  lew  Reynolds  number,  the  empirical 
correlation   was    in    the   form: 

Nu    =    Nug  (1+C  (Ear)  m/Gz)  n,  (2.17) 

where   the   values    of    the    constants    were: 

C    =    0. 03  ra    =    0.75  n    =    1/3 

and    where      Nu      =   1 .  86  (RePrd/L)  °«  3 3      for      large  L/d   values, 
and    Nu      =    2.67    Gz0*33      for    small   L/d    values. 

This  was  valid  within  +10^  to  -17%  for  Reynolds  numbers 
between  162  and  2,700  and  rotational  Rayleigh  numbers  to 
20,000,000.  It    was      also      found   that      the      length   of      the 

entrance  region  seemed  to  have  no  affect  on  the  correlation, 
which  would  imply  that  a  secondary  flow  existed  in  both 
configurations.  An  interesting  point  in  this  paper  is  that 
the  rotational  Grashcf  number  (hence  the  rotational  Rayleigh 
number)  was  calculated  using  the  temperature  difference 
between  the  wall  and  the  fluid  instead  of  the  axial  tempera- 
ture   gradient. 

Woods  and  Morris  investigated  laminar  flow  in  the 
rotor  windings  of  directly-cooled,  electrical  machines  in 
197U      [Ref.    38]    and      noted    that      if      water    was      used    as      the 
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coolant  in  lieu  of  hydrogen  (which  is  currently  widely 
used),  an  efficiency  improvement  of  >0.5a  could  be  obtained. 
Given  that  the  0.5~  cf  a  large,  say  1,000  Ktf,  generator  was 
the  only  consideration  made  (and  neglecting  the  size  and 
weight  reduction  and  the  extended  operating  life  possible 
due  to  the  cooler  operation),  they  proceeded  to  analyze  the 
fundamentals  of      the    problem.  They    surveyed      the   previous 

works    directly    applicable    to   this    area    of   cooling; 

1)  Morris  [  Ref .  14]  used  a  series-expansion  technique 
(valid   fcr   low    rotation   rate   and   low    heating   rates) ,    and 

2)  Mori  and  Nakayama  [ Bef .  28]  assumed  a  secondary  flow 
(which  was  claimed  to  be  valid  fcr  high  rotation  rates) ,  and 
used    an    integral-type  analysis. 

Woods  and  Morris  claimed  both  were  inadequate  due  to  the 
restrictions  imposed  fcy  the  nature  of  their  solutions.  They 
attempted  a  versatile  and  "exact"  solution  for  laminar  flow 
in  the  fully-developed  region  by  solving  the  governing  equa- 
tions wi^h  a  numerical  procedure.  Numerical  solutions  were 
presented  for  both  the  friction  factors  and  the  Nusselt 
numbers.  Their  analytical  values  were  compared  with  experi- 
mental results  and  to  correlations  of  previous  works. 
Eiscr epancies  were  explained  and  the  difficulty  in  obtaining 
fully-de velcped  conditions  in  the  experiments  of  their  work 
and   previous      works      was   emphasized.  The      axial      density- 

variaticn  effects  were  noted,  also,  as  a  potential  for  error 
in  their  analysis  at  high  values  of  rotational  P.ayleigh 
numbers. 

Nakayama  and  Fuzicka  evaluated  the  generator 
problem,  where  reasonably  large  radii  of  rotors  (approxi- 
mately 1  m)  and  high  RPM  (3600)  create  high  centrifugal 
acceleration  affects  [Ref.  39].  This  acceleration  causes 
the    centrifugal    buoyancy      term   to   be    significant.  It   also 
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emphasizes  the  Coriolis  term  in  the  radial  tube  segments  and 
curved  inlet  regions.  They  presented  correlations  for  the 
ratio  of  friction  factors  and  the  ratio  of  Nusselt  numbers, 
arriving  at  friction-factor  ra-ios  for  each  of  the  three 
regions  (radial,  curved  inlet,  and  axial)  ,  and  compared  the 
results  to  experimental  data.  Their  correlations  are  listed 
telo v  : 

f/f0    =    2.2Ro-o.33       (radial)  (2.18) 


f/f0     =1.5  Ro-o-3       (ceil    end)  (2.19) 


f/fo     =    14.0    Ro-o- 8       (axial)  (2.20) 


Nu    =    Re0-3   Pro- ♦  {0.01U (Re/Ro2-5) o -12*}        (radial)        (2.21) 

Their  report  confirms  the  results  of  Nakayama's  correlation 
(equation    (2.14))    for  the   axial    sections. 

The  application  of  the  technology  of  rotating 
cooling  schemes  is  net  limited  to  the  printing  industry  and 
electric  machine  area;  extensive  research  has  been  devoted 
to  this  field  of  study  by  the  gas  turbine  industry.  The 
Advisory  Group  for  Aerospace  Research  and  Development 
(AGARE)  held  a  meeting  in  September,  1977,  at  Ankara, 
Turkey,  to  discuss  High  Temperature  Problems  in  Gas  Turbine 
Engines.  Numerous    papers      were   presented      on    the      subject 

including  a  paper  by  W.D.  Morris  on  flow  and  heat  transfer 
in  rotating  coolant  channels  [Bef.  40  ].  He  used  a  selection 
of      exDerimental      results      to    illustrate      the      influence     of 
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rotation  on  heat  transfer,  and  demonstrated  that  Coriolis 
a:\l  centripetal  inertial  effects  significantly  alter  the 
heat -transfer  characteristics  compared  to  the  stationary 
case.  Using  the  governing  equations,  he  illustrated  how  the 
Navier-Stokes  equations  are  modified  to  include  the  acceler- 
ation effects  of  the  Coriolis  forces  and  the  centripetal 
forces.  For  constant  properties,  these  equations  than  show 
how  at  distances  sufficiently  downstream  from  the  entry 
region,  where  the  axial  gradients  of  the  velocity  are  negli- 
gible, the  elimination  of  the  pressure-gradient  terms  from 
the  radial  and  tangential  momentum  equations  causes  the 
Coriolis  terms  to  vanish  identically  as  a  source  for  the 
creation   of   secondary  cross-stream   flow;    see   also   [Ref.    14]. 

In  the  entry  region,  the  Coriolis  terms  interact 
with  the  developing  axial  velocity  and  cr=ate  secondary 
flows  perpendicular  to  the  axial  direction,  even  with  the 
constant    property      condition.  Finally,  in      the      fully- 

developed  region,  the  effect  of  the  rotation  manifests 
itself  only  as  a  cross-stream  pressure  distribution  if  buoy- 
ancy terms  are  neglected.  If  the  density  is  allowed  to  vary 
with  temperature,  the  centripetal  acceleration  terms  need 
only  be  included  as  the  Coriolis  terms  do  not  affect  this 
action.  This  cross-stream,  buoyancy-induced  flow  gives  rise 
to  greatly-enhanced  heat  transfer  and  an  attendant  increase 
in  flow  resistance.  The  effects  of  the  eccentricity,  with 
regard  tc  the  Coriolis  terms,  has  little  effect,  provided 
that  the  ratio  of  the  radius  (H)  to  the  diameter  of  the  tube 
(d)  is  greater  than  5  (only  very  small  radius  of  rotation 
will    have   any   noticeable  effects). 

Harto  [Ref.  41]  lists  the  papers  presented  at  the 
14th  Symposium  of  the  International  Centre  for  Heat  and  Mass 
Transfer  (ICHHT)  held  in  Dubrovnik,  Yugoslavia,  in  1982. 
Papers  of  particular  interest  to  this  study  included  a  paper 
by   Jchnscn      and    Morris     [Ref.    42],      wherein      the    concept     of 
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centripetal-acceleration-induced  secondary  flow  was  proves 
to  be  cf  little  importance,  and  its  influence  was  show?,  to 
be      totally  hydrostatic.  Their      work      also    confirmed      the 

effect  cf  Coriolis  acceleration  with  regard  tc  secondary 
flows,  especially  in  sections  with  relatively  small  L/d 
ratios.      v7o   quantitative  recommendations   were   presented. 

5 .      Flow  Transit  ions 

Mori  [Ref.  25]  states,  when  considering  turbulent 
flow,  that  it  is  not  necessary  to  consider  the  influence  of 
buoyancy  on  heat  transfer.  The  secondary  flow  was  found  to 
suppress  the  turbulence  level  when  the  turbulence  at  the 
inlet  region  was  high  and  an  empirical  formula  for  the  crit- 
ical Reynolds  number,  in  terms  cf  the  Reynolds  number  and 
Rayleigh   number    product,    was   calculated: 

Re         =    123(ReRa)  o-  2s.  (2.  22) 

cr 

Contrary  to  this,  when  the  turbulence  at  the  inlet  region 
was  low,  the  critical  Reynolds  number  was  higher  (7,700  vs. 
2,000),  and  heating  decreased  the  critcal  value.  The  net 
result  was  that  when  the  ReRa  value  was  high,  the  secondary 
flow  caused  by  buoyancy  makes  the  critical  value  of  Reynolds 
number  tend  toward  the  same  value,  whether  the  turbulence  in 
the  inlet  is  high  or  low.  In  the  case  of  low  turbulence  in 
the  inlet,  the  value  for  the  critical  Reynolds  number  was 
given    by: 

Re        =   Re        /(HO.  14  ReRa*10~s).  (2.23) 

cr  crg 

The  graph  shown  as  Figure  62  in  the  Krieth  article  [Ref.  13] 
is  shewn  below  as  Figure  2.1.  This  graph  shows  the  heat- 
transfer  as  a  functicn  of  both  the  RPM  and  the  flow-through 
Reynolds   number.         It   is   easy    to    see    the    flow's   transitional 
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deoendence    en   the   Reynolds    number   as    a   function    of    PPM.      For 


:ha    experimental   model  being 
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Rar=*7,40C,      -he   transition      is   delayed   to   Re=11,000      Eor    the 
ccrrelation    that    assumes    high    turbulence   at   the    entrance. 


Figure   2.1.      Effect   of    Rotational   Speed 

*  in  a   Rotating   Pipe. 


The  transition  to  turbulent  flow  is  suppressed  by 
rotation  and  has  been  qualitatively  noted  by  numerous 
authors.  An  exact  correlation  for  transition  to  turbulent 
flow  in  a  r eta  ting-reference  frame  with  heating  is  net  known 
to  exist.  From  the  available  literature,  it  is  clear  that 
the   transition    follows  a   path      from    the   laminar   correlations 
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to  the  turbulent  correlations.  This  path,  in  the  case  of 
the  rotating  reference,  has  a  "dip"  as  shewn  qualitatively 
ir.  a  plot  of  the  Nusselt  number  against  the  Reynolds  number. 
Figure  2.2,  whereas  in  the  non-rotating  cass,  no  "dip" 
occurs. 
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6 .      Summary 

Eased  upon  the  information  reviewed  in  the  litera- 
ture, it  may  be  concluded  that  the  three  regions  within  the 
motor/generator  application  (shown  in  Figure  1.11)  hava  very 
different   heat-transfer   characteristics. 

1)  In  the  radial  section  of  the  coolant  path,  the  Coriolis 
forces  tend  to  dominate  the  heat  transfer  and  friction. 
This  effect  is  to  enhance  the  heat  transfer  and  increase  the 
pressure  drop  in  these  sections.  The  heat-transfer  enhance- 
ment is  substantial  and  could  be  effectively  exploited  for 
the  heat  rejection  from  the  device  by  axial  fcrced-air 
convection    of    the    substantially-hollow,    synchronous    rotor. 
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2)  In  the  coil-end  regions  (both  at  the  radial  connection 
points  and  the  coil- end  loops)  the  Coriolis  force  is  strong 
and  is  coupled  with  the  strong,  centripetal  forces.  The 
result  is  enhanced  heat-transfer  and  fricticr.-f  actor 
effects,    but  to    a   lesser   degree  than   in   the   radial    sections. 

3)  Irs  the  winding  bars,  the  centripetal  forces  dominate: 
these  forces  result  in  hydrostatic  effects  which  limit  flows 
that  would  be  affected  by  the  Coriolis  forces  and  result  in 
some  degree  of  secondary  flew  due  to  buoyancy.  The  amount 
cf  secondary  flow  is  questionable  and  is  somewhat  affected 
by  the  degree  of  turbulence  and  the  amount  enhancement  is 
less    specific. 

The  enhancement  in  the  laminar-flow  heat  transfer 
has  teen  reported  by  numerous  authors.  This  can  be  seen  in 
Figure        2.3.  Both        ths      Nakayama        [Ref.    29]      and        the 

Woods-Morris  [Ref.  31]  correlations  for  Nusselt  number  fall 
between  the  classical  laminar  value  of  48/11  and  the  turbu- 
lent value  cf  the  Dittus-3o elter  correlation.  They  have  the 
same  trend  throughout  the  range  of  RPM  shown  and  differ  by 
only  a  few  percentage  points.  Apparently,  the  effect  of 
rotation  in  laminar  flow  enhances  the  heat  transfer  because 
cf  intense  secondary  flows  due  to  centripetal  forces.  An 
interesting  point  is  that,  at  the  higher  RPM  values,  the 
correlations  of  Nakayama,  Woods-Morris,  Stephenson,  and 
Dittus-Bcetler    all   are   exceptionally   close    together. 

The  heat -transfer  predictions  for  turbulent,  fully- 
develcped,  parallel  flows  are  shown  in  Figure  2.4.  This 
figure  illustrat es  that  the  resultant  heat  transfer  is  some- 
where in  the  vicinity  of  the  classical  turbulent  correlation 
cf  Dittus-Bcelter   for  the   current   model. 
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The  hea*--  transfer  is  either  somewhat  enhanced  (on  the  order 
cf  30*)  in  the  lower  RPM  ranges  (Nakayama  (  Ref .  29]  )  ,  or 
somewhat  suppressed  by  the  high  centripetal  acceleration 
forces  as  shown  by  the  correlation  of  Stephenson  [Ref,  43]. 
The  enhancement  of  the  heat  transfer  in  the  Nakayama 
correlation  [Ref.  29]  in  the  lower  RPM  range  is  assumed  to 
he  due  to  the  residual  Coriolis  effects  from  the  entrance 
regions  remaining  longer  in  the  lower  centripetal  accelera- 
tion   field. 

B.       HISTORICM.    DEVELCPHENT    OF    THERMOSYPHON    COOLING 

The  thermosyphon  is  a  device  that  creates  the  motion  of 
the      working  fluid      solely    by      buoyancy   forces.  It   has     a 

region  wher<=  heat  is  supplied  to  the  working  fluid  and  a 
region  where  heat  is  rejected.  The  less  buoyant,  cooler 
porticn  cf  tha  fluid  moves  in  the  direction  of  the  accelera- 
tion vector,  while  the  opposite  is  occurring  for  the  warmer 
portion.  Th=  thermcsyphon  also  has  a  third  region,  the 
coupling  region,  with  problems  of  fluid  shear,  mixing  and 
entrainment.  Figure  2.5  illustrates  the  conceptual  arrange- 
ments  of   the  three    basic  thermosyphon    systems. 

Kreith  [Ref.  13]  mentioned  in  his  review  that  the  ther- 
nosyphon  (i.=.,  a  liguid -filled  container)  had  a  great 
potential  for  cooling  of  slectric  rotating  machinery.  In 
motors  and  generators,  the  use  of  a  thermosyphon  would 
preclude  the  necessity  of  using  rotating  seals,  which  have 
marginal  reliability  in  applications  of  high  RPM  and  volu- 
metric flow.  In  1973,  Japikse  [Ref.  44]  gave  a  thorough 
review  of  the  literature  regarding  thermosyphon  technology. 
Beginning  with  a  reference  to  the  Davies  and  Morris  paper 
[Ref.  45  ]f  which  will  be  discussed  in  detail  later,  Japikse 
categorizes   ther mosy phons  as   follows: 
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1)  the  nature  of  boundaries  (either  opened  or  closed  to 
mass    flow) , 

2)  the  regime  of  heat  transfer  (purely  natural  convec- 
tion   cr    sixed   natural  and   forced   convection), 

3)  the   number   and  types    of   phases   present,   and 

U)  the  r.ature  of  the  body  fores  present  (gravitational 
cr   rotational)  . 

These  broad  categories  are  than  restriced  to  M....a 
prescribed  circulating  fluid  system  driven  by  thermal  buoy- 
ancy   forces." 

The  use  of  the  closed  thermosyphon  in  rotating  refer- 
ences, such  as  in  gas  turbine  cooling,  has  been  studied  in 
detail.  In  1965,  Eayley  and  Lock  [Ref.  46]  experimentally 
investigated  the  performance  of  closed  thermosyphons.  They 
utilized  the  theoretical  development  of  the  open  thermosy- 
phon, modified  the  theory,  and  experimentally  verified  their 
results.  It  was  shewn  that  the  critical  operating  parame- 
ters   include: 

1)  Length -to- Diameter  Ratio.  This  ratio  controls  the 
characteristics  of  th€  flow  and  determines  the  extent  of  the 
coupling   region. 

2)  Heated-L*ngth-to-Cooled-Length  Ratio.  As  this  ratio 
tends  toward  zero,  the  analysis  is  identical  to  the  open 
ther mesyphon  correlations.  This  restricts  the  use  of  the 
analysis  for  the  clcsed  thermosyphon  to  Lh/Lc  ratio  to  be 
not    much   greater   than  one. 

3)  Coupling  Region.  This  is  the  region  that  compli- 
cates the  analysis  and  causes  the  Prandtl  number  to  affect 
the  heat  transfer.  It  exists  in  three  distinct  modes; 
conduction,    convection,      and    mixing.         These   are    regarded   as 
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ideals  and  actual  coupling  usually  consists  of  a  combination 
of   ail   three, 

Thermosyphon  technology  was  also  applied  to  transformer 
cooling,  nuclear-reactor  cooling,  heat-exchanger  fins,  home 
and  industrial  furnaces,  cryogenic  cool-down  apparatus, 
steam  tubes  for  bakers*  ovens,  internal  combustion  engine 
cooling,  and  environmental  control  in  the  space  vehicle 
program,  as  well  as  many  others  mentioned  in  the  Japikse 
article. 

In  1971,  Bayley  and  Martin  [Ref.  47]  also  reviewed  the 
state-cf-the-art  technology,  with  particular  emphasis  on 
gas-turbine   applications.  They   studied    both  the      open   and 

closed  thermosyphon  systems.  The  use  of  an  open  thermosy- 
phon system  in  a  rotating  reference,  such  as  in  gas  turbine 
cooling,  as  well  as  in  the  cooling  of  electric,  rotating- 
machir<=ry  area,  has  problems  regarding  the  fluid  selection 
and  its  saturation  temperature  and  pressure.  In  the  open, 
rotating  reference  system,  the  region  that  receives  the  heat 
is  under  high  acceleration  forces  with  high  resultant  pres- 
sures. As  the  natural  flow  due  to  buoyancy  occurs,  the 
heated  fluid  moves  from  a  high  pressure  to  a  relatively  low 
pressure  area  where  a  phase  change  to  vapor  may  occur, 
possibly  even  explosively,  blocking  the  flow.  They  intro- 
duced the  idea  of  using  two-phase  cooling  within  the  closed 
thermosyphon  as  a  way  of  exploiting  this  phenomenon,  by 
using  the  high  heat  fluxes  associated  with  the  latent  heats 
of  evaporation  and  condensaticn  and  the  lower  temperature 
gradients  associated  with  the  phase  change  processes.  This 
also  has  an  added  benefit  of  reduced  weight  over  the 
single-phase     liquid     system.  They        noted      that      further 

research  in  the  area  of  the  Coriolis  force  effects  on  the 
heat -transfer  chacteristics  in  a  rotating  reference  needed 
to    be   accomplished. 
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applying  the  technology  of  the  t her mo syphon  to  electric 
■achines,  Morris  and  Davies  [Hef.  45]  adapted  the  principles 
of  the  closed  thermcsypbon  to  a  clcssd  rotating  loop  that 
could  be  placed  into  the  rotor  assembly  of  an  electrical 
device.  This  configuration  avoids  the  entry  chokina  and 
nixing  that  occurs  in  the  typical  closed  ther mosyphon  and 
generally  avoids  the  complex  mid-tube  exchange  process,  as 
well  as  the  adverse  core-boundary  layer  interaction.  The 
strong  acceleration  force  present  in  the  periphery  of  the 
rotor  of  electric  rotating  machinery,  coupled  with  the 
heating  due  to  the  electrical  losses,  induces  stong  loop 
circulation.  This  beat  can  then  be  rejected  either  to  the 
shaft  of  the  rotor  acting  as  a  heat  exchanger  itself,  tc  the 
air  passing  through  the  hollow  rotor  region  and  the  air 
passing  through  the  end-bell  region  of  the  roxor  ends,  cr  to 
both.  The  shaft,  acting  as  a  heat-exchanger  itself,  can  be 
cooled  with  an  internal,  rotating  heat-pipe  (or  two-phase 
ther mosyphon)  which  could  transfer  the  heat,  axially  to 
another  heat  sink  outside  of  the  motor/generator  casing.  As 
an  alternative,  the  shaft  cculd  be  internally  cooled  by 
forced  convection  of  water  (cr  some  other  fluid)  to  an 
external  heat  exchanger  through  the  non-coupling  end  of  the 
shaft  by  the  use  of  an  axial  rotating  seal.  The  reliability 
and  the  casualty-control  aspects  of  the  axial  rotating  seal 
make  it  preferable  tc  the  radial  rotating  seal  mentioned 
earlier. 

Theoretical  analysis  of  the  loop-flow  has  used  simple, 
or.e-dimensicr.al  force  and  energy  balances,  and  need  net  be 
repeated  here.  The  closed  loop  is  shown  in  Figure  2.6  as 
conceptually  applied  to  a  rotor.  Of  course,  numerous  loops 
would  be  designed  into  the  rotor  to  carry  away  the  heat 
generated.  Incorporation  of  a  one-way  valve  was  suggested 
as  being  necessary  tc  insure  one-directional  flow  created  by 
the    mean   temperature    difference    and   the    pressure    differences 
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Figure  2.6    Rotating  Loop  Theraosyphon  for  a  Rotor. 

in  the  limbs  of  the  circuit.  Experimental  testing  was 
conducted  in  a  test  rig  similar  to  the  rig  proposed  for  the 
evaluation  of  the  systems  studied  herein  and  now  under 
construction  at  DTNSBDC.  A  modest  range  of  rotation  rates 
was  employed  in  the  Wcrris  and  Davies  test  rig,  typical  of 
most  meter  applications.  The  details  of  their  tes-1-  section 
were  very  similar  tc  the  model  used  for  analysis  herein: 
25.1  mm  (one  inch)  diameter  copper  rod  through  which  a  6.35 
mm  (0.25  inch)  hole  was  bored.  The  experimental  results 
were  placed  in  the  form: 
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Nu    =   f  (Grr,Ro,Pr,AC).  (2.2^ 

The  ratic  of  angular  acceleration  to  gravitational  accelera- 
tion (Ac)  was  the  authors1  way  cf  including  the  rotational 
effects  (previous  authors  have  used  the  rotational  Grashof 
number).      The   final    correlation   was: 

Nu    =    0.1505    Aco.735   pa2.*s    pr/Gr  (2.25) 

for   the    range   of    103<  NuGr/Pr   <106. 

The  critical  consideration  in  the  analysis  of  the 
closed-loop  t hermosyphon  is  the  temperature  rise  in  the 
heated    limb    (conductor  bar)  .      The   model    described    in    Chapter 

I  was  analyzed  by  using  the  governing  relationships  from  the 
Davies  and  Morris  [Ref.  45]  paper  iterated  against  the  Woods 
and  Morris  [Ref.  31]  correlation  for  laminar,  fully- 
develcped  flow  in  a  heated  circular  tube  rotating  about  a 
parallel  axis.  It  can  be  seen  from  Figure  2.7  that  the  heat 
transfer      increases        with      RPM       (the        higher     acceleration 

TABLE    II 
Temperature  Rise  in  the   Closed-Loop   Thermosyphon 

RP M  J.em£.    Rise  (   C)  * 

300  209.0 

1000  93.5 

3000  41.0 

3600  35.5 

5000  27.0 
*    for    a   cold-side    temperature   of   45    C. 

effect).  For  the  conditions  on  the  model  which  constructed 
Figure    2.7,    the    values  of   temperature    rise    is   shown    in    Table 

II  for    various    RPM's.        This      clearly    shows   that    this    method 
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cf  cooling,  if  selected,  would  be  inappropriate  for  the 
■otoi  device  (i.e.,  lew  RPH) ;  on  the  other  hand,  it  could 
be  very  reasonable  for  the  generator.  The  friction  factor 
was  assumed  to  be  6  4/Re  and  no  correction  was  mads  for  the 
rotation  in  the  friction-factor  calculation;  the  minor 
losses  were  neglected.  If  the  friction  factor  is  increased 
due  tc  the  rotation,  as  in  equation  (2.U),  and  minor  losses 
(laminar)  are  included,  the  temperature  rise  will  be 
increased   slightly.  The   trend      will   be      the   same      and   the 

temperature  rise  at  the  higher  RPM's  will  still  be  accep- 
table. The  actual  design  of  the  closed-loop  thermosyphon 
will  require  an  optimization  of  the  heat-transfer  calcula- 
tions and  the  friction- fact  or  calculations.  This  optimiza- 
tion may  indicate  that  a  larger,  single,  ciosed-t her mosyphon 
loop  for  a  conductor  bar  group  is  preferable  over  individual 
conductor   loops. 

C-       ADVANTAGES/DISADVANTAGES    OF    LIQUID    COOLING 

The  main  objectives  cf  cooling  are  to  remove  the  heat 
created  by  the  electrical  losses  and  improve  the  operating 
efficiency,  expend  the  lifetime  of  the  insulation,  and 
reduce      the   overall      weight    and      size   of      the   device.  The 

feature  cf  liquid  coding  (vs.  gas  cooling)  is  the  superior 
capacity  of  liquids,  especially  water,  to  remove  heat  with  a 
small  rise  in  temperature.  The  principle  disadvantage  of 
using  liquids,  or  gaseous  fluids  other  than  air,  is  the 
desiqn  cf  the  piping,  coupling,  and  secondary  heat  exchanger 
required  to  support  the  device  cooled  by  the  fluid.  Since 
one  of  the  main  objectives  is  to  reduce  weight  and  size,  the 
fluid  chosen  must  be  able  to  remove  enough  heat  per  unit 
mass  to  accomplish  this  goal  within  the  device. 
Additionally,  since  the  application  under  evaluation  is 
marine    combatants,    the   fluid   should   be   non-toxic,    and    thc 
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systea  should  b?  durable  enough  to  withstand  damag?  and 
maintain  operational  safely.  The  popular  fluid  tc  date  in 
land-based  installations  has  been  hydrogen  which  is  clearly 
an  unacceptable  choice  for  the  marine  combatant  due  tc  its 
explosive  nature  and  since  the  size  reduction  obtained  is 
net      sufficient      to      warrant   system      development.  Of      the 

remaining  fluids,  the  obvious  choice  for  its  high  heat 
capacity,  availability,  and  safety,  is  water.  The  use  of 
water  has  been  tested  both  in  staters  and,  on  a  limited 
scale,    in   rotors. 

The  principle  problem  to  overcome  in  the  use  of  water  is 
in  the  rotary  seals  needed  to  direct  the  fluid  into/cut  of 
the  machine.  Additionally,  when  the  system  inside  the 
device  is  analyzed,  the  ability  of  a  design  to  withstand 
damage      is    questionable.  The   design      would    be      critically 

dependent  on  the  fluid  within  the  rotor  for  operation,  even 
in  a  casualty  mode.  The  supply  of  cooling  water  to  a  stator 
would  net  be  as  critical,  since  any  casualty  could  be 
repaired  easily  with  shipboard  damage-control  equipment,  and 
emergency  cooling  of  a  stator  could  b«  accomplished  by 
external  means  due  to  the  thermal  conductivity  of  the  stator 
housing.  For  the  rctor,  this  leaves  a  choice;  either  not 
utilize  the  enhanced  cooling  and  suffer  the  weight  and  size 
penalties,  or  develop  a  closed  system  for  cooling  the  rotor. 
The  closed  system  for  the  rotor  is  not  necessarily  isolated. 
A  closed-lccp,  thermesyphon  could  be  placed  in  the  rotor, 
transferring  the  heat  from  the  rotor  to  the  rotor  shaft. 
This  loop  could  have  additional  heat  transfer  to  air  passing 
through  the  core  of  the  rotor.  Further,  the  heat  trans- 
ferred tc  the  shaft  could  be  removed  to  the  ambient  air  or 
an  external  heat  exchanger  utilizing  the  non-coupling  end  of 
the  motor  by  the  use  of  circulating  water  through  the  shaft, 
or   by   a    rotating   heat  pipe    as    seen   in   [Hef.    4  8]. 
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An  additional  method  fcr  cooling  the  rotor,  which  has 
beer,  studied  by  numerous  authors,  is  the  use  of  heat  si  pes 
cr   closed,    two-phase   thermo syphons.        In   this   case,    the   heat 

is  transferred  from  the  rotor  to  air  circulated  through  the 
end-bell  areas  of  the  device.  This  method  of  rotor  cooling 
greatly  enhances  tie  survivability  of  the  device  in  a 
damage-control  sense.  Each  heat  pipe  is  independent  of  the 
others  and  damage  to  a  few,  unless  all  in  the  same  group  of 
adjacent  conductor  bars,  would  not  jepardize  the  device; 
indeed,  even  if  all  were  in  the  same  group,  limited  opera- 
tion   could   continue. 
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III.     TWO-PH aSS    FLOID    COOLING 

The  concept  of  utilizing  the  latent  heat  of  vaporization 
to  transfer  heat  is  widely  used  in  heat  exchangers  every- 
where, from  household  and  industrial  boilers  and  space 
heaters  to  marine  propulsion  units.  Air-conditioning  equip- 
ment and  heat-pump  units  also  employ  this  idea  to  advanta- 
geously utilize  the  low  temperature  of  vaporization  of 
fluorc-carbcn  compounds  to  remove  the  heat  from  the  ambient 
and  release  this  through  condensation  under  pressure,  thus 
transferring  the  heat  in  the  desired  direction  at  the  cost 
of  the  energy  to  pump  and  pressurize  the  gas  to  the  point  of 
condensation.  In  each  of  these  two-phase  applications,  the 
transport  of  the  working  fluid  may  require  a  pump,  or  it  may 
be  accomplished  by  the  buoyant  force  that  exists  in  an 
acceleration  field,  such  as  that  "due  to  gravity.  This  is 
the  principle  that  is  employed  in  the  closed  two-phase  ther- 
mosyphon.  This  device  is  similar  to  the  classical  "heat 
pipe"  and  the  only  difference  is  that  the  "pumping"  method 
in  the  heat  pipe  uses  capillary  action  instead  of  accelera- 
tion forces  to  pump  the  fluid.  Thus,  a  "heat  pipe"  used  in 
a  rotating  reference  utilizing  the  acceleration  forces  to 
transport  the  liquid  from  the  condenser  section  to  the 
evaporator  section  (in  lieu  of  a  wick)  is,  technically,  a 
two-phase  ther mosy phcn.  Figure  3.1  illustrates  the  opera- 
tion of  the  closed  two-phase  thermosyphon .  The  evaporator 
end  must  be  in  the  direction  of  the  acceleration  field  for 
the  liquid  to  be  transported  from  the  condenser.  Figure  3.2 
illustrates  the  operation  of  the  heat  pipe,  which  utilizes 
the  capillary  action  of  the  wicking  material  to  transport 
the      liquid    from      the     condenser      to    the      evaporator.  The 

remainder    of  the      chapter  shall   use    these      definitions    for  a 
discussion    of   both    devices. 
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Figure    3.1         Typical   Two-Phase   Thermosyphon. 
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Figure    3.2        Typical   Heat   Pipe. 

A.       HISTORICAL    DEVELCPHEHT 

The  development  cf  the  heat  pipe  (which  will  apply 
equally  to  the  development  required  for  the  two-phase  ther- 
mosyphon) was  begun  by  Gaugler  [Ref.  49  ]  in  1944  and  by 
Trefethen      [Ref.    50]    in      1962.  Trefethen      was    working     on 

cooling      methods   for      spacecraft       (zero-g   environment)         for 
General      Electric     Company      and      discovered      that      capillary 
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pumping  can  be  a  very  valuable  area  for  further  development. 
Working  independently,  Grover,  ex  al.  (Ref.  51]  reinvented 
this  concept;  he  widely  published  the  first  applications  of 
the  device  and  gave  it  i~s  name.  Kraus,  et  al-  [Rex.  52] 
have  reported  this  historical  development  and  have  given 
sample  calculations  for  the  design  of  a  heat  pipe.  Chi 
[Ref.  53]  also  reviews  the  theoretical  development  of  heat 
pipes,  and  provides  the  design  procedures,  including  sample 
calculations,    for  their   utilization. 


B.       CONVENTIONAL,    CLOSED   TWO-PHASE    THERMOSYPHON    ANALYSIS 

As  previously  mentioned,  the  appealing  factors  associ- 
ated with  the  two-phase  utilization  include:  high  heat 
fluxes  associated  with  phase  changes,  lower  temperature 
gradients  associated  with  these  processes,  and  reduced 
weight  of  the  two-phase  system  over  the  single-phase  liquid 
system. 

Research  by  Cohen  and  Bayley  [Ref.  54]  (referred  to  and 
discussed  by  Japikse  [Ref.  44  ]  )  found  that  the  amount  of 
liquid  filling  the  system  functionally  affected  to  the  heat 
transfer  in  both  rotating  and  static  tests.  They  found  that 
heat  transfer  increased  as  the  percent  liquid  in  the  system 
was  increased  to  approximately  1.5%  by  volume,  then 
decreased  to  an  intermediate  value,  and  finally  increased  to 
approximately  the  same  value  as  in  the  1.5%  case.  This  was 
related  to  tha  following  process:  in  the  low  filling  situ- 
ation, condensate  returns  to  the  evaporator  section  and 
forms  a  thin  film  on  the  walls.  It  is  in  this  film  that  the 
heat  transfer  occurs  and  the  phase  change  takes  place.  In 
the  case  of  the  co mp letely- filled  evaporator,  a  liquid  pool 
region  develops,  and  a  condition  of  nucleate  boiling  exists. 
If  the  filling  is  insufficient,  with  regard  to  the  heat  flux 
being   transported,       the      pool/film    will   not   continue      to  the 
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end  of  the  heated  section  and  "dry -oat"  will  occur.  This 
"dry-cat"  will  assume  a  drop-wise  instead  of  film-wise  phase 
change.  The  geometry  of  the  regions  in  which  the  films  are 
likely  tc   form,     influences    the   results. 

Lee  and  Hital  ( Eef .  55]  conducted  experiments  with  an 
electrically-heated,  water-cooled  thermosyphon  using  water 
and  freon  as  the  test  fluids  and  varied  the  filling  guan- 
tityr  Lh/Lc,  pressure  (or  I  ),  and  heat  flux  (L  /L  is  the 
ratio  cf  the  length  of  the  heated  section  to  the  cooled 
section) .  The  result      of      the      filling    guantity      on      heat 

transfer  was  the  same  as  that  reported  by  Cohen  and  Bayley: 
increasing  heat  transfer  with  filling  to  a  point  and  then 
decreasing  beyond  that  value  to  an  intermediate  value  and 
increasing  to  the  case  of  the  completely-filled  evaporator. 
The  effect  of  decreasing  1^/1  was  to  increase  the  heat 
transfer  within  the  range  0.8  to  2.0  and  the  advantage  of 
larger  ccndenser  area  was  evident.  The  heat -transfer  coef- 
ficient was  found  to  increase  significantly  with  increasing 
mean    pressure   due   tc   at   least   three    factors: 

1)  since  the   mass    flew    for    a    given   heat   flux  is    nearly 
constant  and   density  of    the   liguid   increases    with 
pressure,    a    lower   pressure   drop    (and   lower     AT)    is 
necessary   for   the   same    heat    flux; 

2)  for    a   larger    PM4.  ,    P   _+.    varies    much    more   rapidly   with 

sat         sat 

T  for  the   fluids   considered,    hence   reguiring   smaller 

sat 
AT*s   at   higher    pressures;    and 

3)  for  lower  pressure  drops,  a  more  favorable  force  balance 
exists  on  the  condensate  film  permitting  a  faster  liguid 
return. 

water  was  found  to  give  heat  fluxes  superior  to  those  of 
freon,  for  for  the  same  aT,  due  to  the  larger  values  of 
latent   heat   of    vaporation,    and   thermal   conductivity. 
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Lee   and     Hital   [  Bef  •    55]  also  considered     the    analytical 

problem  cf  predicting  the  maximum  heat-transfer  rate  for  a 
laminar  film,  constant-wall-temperature  condenser  and  a 
ccnstant-heat-f lux  evaporator.  Neglecting  the  forces  due  to 
vapor  pressure  drop  and  momentum  changes,  and  using  a  force 
balance  on  the  falling  film,  balancing  the  effects  of 
gravity  and  fluid  shear,  they  related  the  mass  flow  rate  to 
the  heat  flux.  Using  a  local  energy  balance  and  an  overall 
energy  balance  between  the  condensing  section  and  the  evapo- 
rator section,  the  following  relations  for  the  heat  trans- 
ferred      (q)      and      the     saturation      temperature    (T        )         were 

sat 
obtained : 


q    =    (p2R3h   3/2u    L,g  }C  (3.1) 

f  g  "   K- 


and 


q    =     {k  (T    -T    )/(PL,   /L    )}C/D  (3.2) 

s      c  n      c 

where  C   =    1/8  -     1/2    y*    ♦    3/8    y*    -    1/2    y*    In    y    , 

D  =  y*ln  y  (1/2  in  y2  -1/2)  +  1/8  y*  *y2  (1/2  In  y  -1/4)  =1/8 
and   y   =    1   -      o/r.  The  simplifying   assumptions    (forces   due 

to   vapcr      shear    and    momentum      changes    were    neglected      in   the 
force    balance)       cause  an  error    (of  as    much   as  a    factor    cf    2) 
to   occur,        but    qualitative      behavior    is      correct    as      far    as 
Lv,/L_  ,      ?!      .and    working   fluid   are   concerned.        This   develop- 

n        C  Sat 

ment  is  also  included  in  the  survey  by  Japikse  [Ref.  44]. 
This  heat  transfer  is  graphically  presented  in  the  Lee  and 
Hital  paper  in  their  Figure  12,  shown  here  as  Figure  3.3. 
This  clearly  shows  that  the  quantity  of  heat  that  can  be 
removed  by  a  two  phase  water  device  is  large,  even  at  a  low 
saturation  temperature.  The  report  also  indicates  that  for 
water  the  increasing  ratio  of  heated  length  to  cooled 
length,  L+,  decreases  the  maximum  heat  transfer.  The  heat 
transfer  is  also   very  sensitive  to  the   operating  pressure; 
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Figure   3.3        Comparison   of    the   Experimental   Results   with   the 
Analytical   Prediction-f rom   Lee   &    Mital   paper. 


the  heat-tr arsf er  coefficient  increases  appreciably  with  the 
mean  operating  pressure.  The  temperature  drop  increases 
with  the  pressure  drcp  alcng  the  length  of  the  t hermosyphon 
tube  (and  the  saturation  pressure  gradient).  As  the  pres- 
sure increases,  the  vapor  specific  volume  decreases, 
resulting  in  a  decrease  in  the  pressure  drcp;  the  mass  flow 
rate  of  the  vapor  is  essentially  constant  for  transferring  a 
given  heat  input  rate.  In  the  previous  papers,  the 
vertically-oriented  thermosyphon  was  considered  and  the 
current  icdel  (for  seme  of  the  orientations)  requires  the 
device  to  be  oriented  in  the  horizontal  direction  (perpendi- 
cular to  the  applied  acceleration  field) .  Although  the  Lee 
and  Mital  paper  was  considering  a  vertically-oriented,  two- 
phase  ther mos7 phon,  the  general  results  are  the  same  for  the 
horizontal    tube. 
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C.       LIHITS    OF    OPERATION 

Hhen  considering  the  3? sign  of  a  two-phase  thermosyphon, 
or  a  heat  pipe,  limitations  to  the  heat  transfer  must  be 
considered;  four  of  these  limits  are  common  to  both  devices 
(sonic,  entrain ment,  toiling,  and  condensing).  A  qualita- 
tive comparison  of  these  heat-tr ansf er  limitations  as  a 
function  of  the  saturation  temperature  is  civ  en  in  fRef«  **8] 
and    is   shown  in    Figur<=   3.4. 


CONDENSING 
LIMIT 


SATURATION     TEMPERATURE 


Figure   3.U        Operating  Limits    of   Rotating   Heat   Pipes. 

The  sonic  limit  and  the  boiling  limit  have  been  readily 
analyzed.  The  entrainment  and  the  condensation  limits  are 
not  as  well-known  and  little,  if  any,  literature  describing 
these  phenomena  exist.  This  is  an  area  in  which  further 
research  is   needed. 

1)         Sonic     limit:    the   vapor      flow   in      a   two-phase      closed 
system    is      limited    to     the    sonic      velocity    at      the    operating 
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pressure    (also      known  as  the      "choking"    limit). 
limit    is   represented   by: 

Q  =    mh        =p      h    ?  h      ,  {^.3; 

max  zg      ■  v      v    s    fg 

where    V      is   the    sonic  velocity  of   the    vapor.        This    limiting 

velocity   may  be    determined      experimentally,      computed    by   the 

relationship     V    s{(3p/3p)    }  l/z,        or        approximated      by      the 

perfect-gas  relationship    {YET        }  */2. 

sat 

2)  Entrain meet  limit;  the  inter facial  shear  between  the 
liquid  and  the  vapor  will  held  the  liquid  (which  is  flowing 
in  the  opposite  direction)  back  and  starve  the  evaporator 
section  cf  liquid.  This  counter-current  flow,  when  the 
relative  velocity  is  large,  causes  the  interface  to  become 
unstable  which  results  in  waves  at  the  interface.  As  the 
vapor  velocity  continues  to  grow,  droplets  of  liquid  are 
formed  at  the  liquid  surface  as  the  shear  force  exceeds  the 
surface-tension  force.  The  formation  of  these  droplets  and 
their  subsequent  entrainment  in  the  vapor  stream  causes  the 
partial      or   total      stoppage    of      the      flow    (dry-out).  This 

phenomenon  is  generally  governed  by  the  Weber  number  (the 
ratio  cf  the  inertial  force  to  the  liquid  surface-tension 
fcrce)  .  The  Froude  number  is  also  used  to  characterize  the 
phenomenon  cf  the  drcp-wise  entrainment  of  the  liquid  in  the 
vapor    stream.  In    the   application    considered      herein,      the 

formation  cf  the  waves  is  considered  unlikely  due  to  the 
extreme     acceleration   field      present.  Since     the    flow      is 

counter-current  and  is  in  the  presence  of  a  high  accelera- 
tion field  the  entrainment  limitation  thus  will  reduce  to 
the      "hold-up"      limit.  It   will      remain      stratified      until 

hold-up        occurs,  resulting  in        evaporator         dry-out. 

Experimental  determination  of  the  exact  correlation  for  the 
hold-up,   or  entrainment,      limit   is   required   for    a    horizontal 
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tube  rotating  around  a  parallel  axis,  thereby  creating  a 
perpendicular  acceleration  field  many  times  greater  than 
cravitv.      The   en trainmert   limit    is    crenerallv    given    hv' 

2max=hfgpvAvVE'  <3-u> 

where  V^  is  when  the  Weber  number  equals  unity,  the  Froude 
number  equals  unity,  or  is  experimentally  found.  In  another 
approach,  Jaster  and  Kosky  [Bef.  56]  used  the  ratio  of  the 
axial  shear  force  of  the  vapor  flow  on  the  liquid  surface  to 
the  gravitational  (accelert ional)  body  force  upon  the  liquid 
(defined  as  "F")  in  order  to  arrive  at  an  appropriate 
"hold-up"  criterion.  They  correlated  experimental  data  as  a 
function  of  this  "F"  value,  as  the  flow  they  were  studying 
transitioned     from    stratified      to   annular.  The   result      of 

their  experiments  was  that  the  flow  was  stratified  for  "F" 
values  below  five  (5)  and  was  fully-annular  above  twenty- 
nine  (29).  The  value  of  "F"  set  equal  to  five  (5)  was 
established  as  the  criterion  for  the  transition  from  strati- 
fied to  annular  for  the  co-current  flow  case.  The  same  type 
of  analogy  was  developed  by  Collier  and  Wallis  [Ref.  57]  who 
balanced  the  inertia  force  and  the  acceleration  force  to 
scale  stratification  effects  where  their  criterion  was  given 
as  j*  (volumetric  flux)  equal  to  0.25.  Both  of  these  refer- 
ences  point   out    the   need   for   further   research   in    this   area. 

3)  Eoiling    liJ5ii:    Again,        as      discussed      in    the      first 

section  of  this  chapter  in  the  general  discussion  on  two- 
phase  thermosy phons,  the  boiling  limitation  in  the  current 
application  is  due  to  the  creation  of  profuse  nucleaticn  at 
the  evapcrator  section  and  a  resultant  vapor  film  on  the 
evaporator  surface  which  insulates  the  evaporator  wall, 
resulting  in  dry-out  and  overheating.  The  boiling  limit  is 
given  by  the  Zu ber-Kutateladze  prediction  as  discussed  in 
[Ref.    48]   and    is  given  as: 
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Q        -0.131   p    »/2hr     {   fi2R(p   -p    )G    }i/S  (3.5) 

max  B  v  fg  1      v 

where    A        is  the  area   in     the   evaporator   section      for    liquid 

B 
heating. 


4)  Condenser  liait;    The      capacity   of      the   condenser     is 

dependent  on  many  parameters,  such  as  the  geometry,  the 
working  fluid,  the  orientation  of  the  vessel  with  respect  to 
the  acceleration  field,  and  the  operating  conditions.  The 
orientation  of  the  vessel,  for  the  model  considered,  is  -hat 
cf  a  horizontal  tube  (perpendicular  to  the  acceleration 
field)  with  internal  condensation.  The  condensing  limit  is 
the  capacity  to  condense  the  vapor  on  the  inside  surface  of 
the  condenser  section,  and  this  is  shown  in  Figure  3.4.  No 
analysis  on  the  outer  surface  is  presented  in  this  thesis. 
Instead,  the  condenser  limit  is  based  on  the  capacity  of  the 
inside  wall  exposed  to  condensation,  using  an  arbitrary,  but 
reasonable,        wall   temperature.  (Generally,      the      largest 

thermal  resistance  occurs  between  the  ambient  and  the 
condenser  cuter  surface.  Therefore,  heat  transfer  is 
limited  by  this  resistance  unless  the  outside  surface  area 
is      substantially      increased      by   the      use     of      fins.)  The 

condensing  limit  is  given  by  an  equation  from  [Ref.  58]  by 
Collier  for  a  horizontal  tube  with  internal  condensation  in 
a  gravitational    field  as: 

Q         =1      ATF    ((p     (P     -o    )gh       k    3)/ 

max      v  c        1      1      v         fg    f 

^I'W^w    "'/'•  (3.6) 

For  the  rotating  reference  being  considered,  the  "g"  term  in 
equation  (3.6)  is  replaced  by  the  centripetal  acceleration 
cf  2  2R .  The  factor  "F  "  in  this  equation,  allows  for  the 
fact    that   the   rate    of  condensation   on    the    statified    layer    of 
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liquid  is  negligable.  The  value  of  "F  "  depends  on  th<= 
ar.gle     i>     and      is    tabulated    in    £Sef.    58].  Its    value   ranges 

from  zerc,  when  the  tuba  is  full  (no  surface  remains  for 
condensation  to  occur),  ro  0.725  ,  when  the  tubs  empty.  The 
angle     $     is   shown   in    rigure    3.5. 


Figure   3.5        Laminar  Condensation    within  a   Horizontal 

Tube    (from   Collier)  . 


1 .      Com  cut at ion    of   Heat-Transfer    Limits 

The  calculation  of  these  limits  is  an  important  step 
in   the    design    process. 

The  calculations  required  to  analyze  the  configura- 
tion of  the  current  model  are  not  all  well  known.  The 
following  appear  to  be  the  most  applicable  to  the  model 
being   evaluated: 

The   sonic  limit    calculation   is      from   basic   principles   as 

shown   by  equation       (3.3)      by  allowing    the      vapor    velocity  to 
increase  to   the    Mach   value    of   unity. 

The      boiling   limit      is      given      by    the      Zuber- Kutateladze 

prediction   as   discussed      in    [Ref.    48]   and   is      given    by    equa- 
tion    (3.5) 


82 


The      condensing      limit      is   given      by      an      equation      from 

[fief.    58]    ry  Collier    in   equation    (3.6). 

The  literature  reveals  no  correlation  yet  for  the  entrap- 
ment limit  in  the  heated,  horizontal  tube  rotating  about  a 
parallel  axis  with  counter- current  flow.  Since  the  majority 
cf   analyses    previously  conducted    for      this    type   of    flow   have 

»—  »  »   v>    W     . .  -  .    —      •>  —  U      — »  ■«.     „     '—      i.    vt  -«»  —    Ulv.     ..  .      t,  w«    —    Ltl  <■—      *  —    o  _-     V      ■-'     -...  krf"    i  .„  -_    —.,       ■  -■_      .,    ._    w     ?S  -  ii 

sonic  and  above  the  condensing  limit,  the  correlations 
previously  mentioned  have  been  plotted  against  the  airour."1-  of 
fluid  that  is  charged  in  the  tube  (percentage  fill)  for  the 
three      reasonably        known      values      (sonic,  boiling,        and 

condensing      limits)       in      Figure      3.7.  Mso   shown      is      the 

entrainment  limit  computed  from  the  Jaster  and  Kosky 
[Ref-    56]  correlation   developed    for   co-current  flow. 

As  was  previously  noted,  the  condensing  limit  is  the 
factor  that  is  the  mcst  limiting  in  the  design  of  the  ther- 
mosyphon.  The  values  of  Qmax  are  increased  to  a  maximum 
based  on  the  length  of  the  condenser,  the  percentage  fill 
and  the  -  T  availatle  between  the  constant  temperature 
(assumed)  wall  and  the  saturation  temperature  of  the  vapor 
condensing  on  it.  As  already  mentioned,  the  resistance  to 
heat-transfer  afforded  by  the  condenser  outer  walls  to  the 
ambient  may  be  quite  significant.  The  ability  to  design  an 
adequate  fin  arrangement  and  supply  cooling  air  in  suffi- 
cient quantities,  while  minimizing  the  windage  losses, 
itself  is  a  substantial  problem.  From  the  data  obtained  in 
this  rough  model,  there  are  some  self-evident  points  to  be 
considered: 

The   length    of  the   condenser,      from   a   heat-transfer    point 

cf    view,    determines    the    heat-removal    potential, 

The    amount    of   liquid   fill    can   be    optimized, 
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The    temperature    at    which   the      walls    of    the   condenser   are 

maintained    and    the    fcemperat ure    at      which    the   ther mosvphcn    is 
operated   are  critical  as   this      temperature    difference    is   the 


driving   force   in   the    hea- 


For  the  model  being  evaluated,  and  neglecting  (as 
the  limiting  value)  fha  entrainment  limit,  the  closed, 
two-phase  ther  mosy  phcn,  operated  at  a  saturation  temperature 
at,  or  above,  10  0°C  and  with  the  condenser  s act ion  held  at  a 
wall  temperature  of  50°  C  with  rotation  at  a  0.381  m  radius 
and  3,600  RPM  can  transfer  the  required  50  W  of  power  with 
the  condensing  limit  being  the  limiting  value.  Figure  3.7 
is  based  on  a  heat  load  of  50  W  to  be  equally  divided 
between  condenser  sections  in  both  end-bell  regions  (i.e., 
approximately  25  W  each  per  condenser)  .  For  the  3600  RPM 
model  the  percentage  fill  may  increase  to  approximately  HQ%, 
and  still  maintain  the  required  heat  flux.  The  effect  that 
the  FPM  has  with  regard  to  the  temperature  rise  in  the 
heated  section,  makes  the  choice  of  the  closed-loop  thermo- 
syphon  poor  for  the  motor  application,  but  quite  suitable 
for  the  generator.  Further  background  on  the  theory  of 
two-phase  flow  is  available  in  publications  by  Wallis 
[Ref.  59]  and  by  Collier  &  Wallis  in  [Ref.  57].  The  feasi- 
bility of  the  design  has  been  shown  by  Corman  and  McLaughlin 
[Ref.  60]  and  by  Groll,  et  al.  [Ref.  61]  through  experi- 
ments with  heat-pipe  cooling  of  electric  motors  in  the 
scheme  primarily  recommended,  although  the  difficulty  of  the 
liquid  transport  at  lew  rotation  rates  was  noted.  A  diagram 
of  the  cooling  scheme  such  as  recommended  and  tested  by 
Cormac  Z  Mclaughlin  and  Groll,  et  al.,  is  shown  in  Figure 
3.6     (reproduced    from   [Hef.    60]  )  . 
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IV.    CONCLUDING    SUMMARY 

As  discussed  in  the  paper  by  Greena,  at  al.  [Ref.  62], 
the  neei  for  electric  propulsion  of  minimal  cost,,  volume, 
and  weight  is  well  known  and  motivated  by  the  inherent 
advantages      of      this      type    of      propulsion.  k      significant 

constraint  to  be  overcome  in  the  development  and  manufacture 
cf  these  systems  is  the  rate  of  heat,  removal  from  these 
devices.  In  1976,  W.  D.  Morris  [Ref.  63]  mentioned  that 
increasing  the  load  per  unit  weight  of  these  devices 
demanded  that  the  inherent  heat  created  by  the  internal 
losses  be  removed  via  improved  cooling  techniques.  icrris 
also  elaborated  on  the  "complex  arrangement"  of  the  ducting 
and  associated  equipment  required  to  accomplish  this 
increased  cooling.  Further,  he  stated  that  "Ther0  is 
currently  insufficient  technical  information  available  with 
which  to  make  confident  predictions  of  the  effect  of 
Coriolis  and  centripetal  acceleration  on  flow  resistance  and 
heat  transfer  in  the  complex  coolant  circuits...".  Since 
then,  Mcrris  and  his  co-workers  have  attempted  to  fill  the 
literature  gap  in  this  area,  and  have  uncovered  a  great  many 
facts  regarding  this  flow  structure;  they  have  also  uncov- 
ered   a   great  many  questions    for   further   research. 

As  shown  as  early  as  1955  in  an  article  by  Sir  Claude 
Gibb  [Ref.  64],  the  principle  factor  in  the  failure  cf  two 
very  expensive  generating  sets  was  the  inability  of  the 
devices   to      dissipate  their      own   heat.  In   the      late   60*s, 

Peach  [Ref.  65]  wrote  of  the  successes  that  the  English 
Electric  Company  had  with  regard  to  liquid  cooling  of  the 
large  generator  rotors  being  used  in  their  500  HW  devices 
and  predicted  that  the  devices  would  reach  1,000  MW  with 
this    technology.         This      was    followed    in    1970      with    a    second 


report  bj  Peach  [Ref.  66]  where  he  outlined  the  growth  q£ 
U.S.  generator  sizes  and  predicted  that  the  liquid-cooling 
technology  would  be  reauired  in  the  near  future  to  support 
the  devices  being  planned  in  this  country  which  exceeded 
1,000  HH.  The  losses  developed  in  the  large-scale,  high- 
power,  superconductive  devices,  currently  being  developed 
for  large-scale  implementation,  require  extensive  use  of 
cooling  at  cyrogenic  temperatures  in  order  to  function  in 
the  super-conductive  state  in  which  they  operate.  A  discus- 
sion of  this  application  is  presented  by  Schwartz  and  Fonsr 
(Ref.  67  ],  in  which  a  number  of  the  schemes  discussed  herein 
are    shown. 

The      previous   chapters      have      discussed      the    ability      of 
external,  f orced-convective        cooling        to        remove        the 

internally-developed  heat  leads  in  electrical,  rotating 
machinery.  The  reduction  in  size,  the  increased  efficiency, 
and  the  prolonged  life  cycle  are  all  very  important  reasons 
to  utilize  advanced-cooling  concepts.  The  use  of  turbulent 
flow  is  feasible  at  all  PPtt's;  therefore,  both  actor  and 
generator  applications  would  benefit  from  forced -•con  vect  ive 
liquid  cooling  utilizing  turbulent  flow.  Laminar-flow  heat 
transfer  is  feasible  only  at  higher  RPM's,  which  tends  to 
limit  the  application  of  this  method  in  regard  to  motors. 
The  actual  transition  from  laminar  to  turbulent  flow  is  not 
well  known  and  further  investigation  is  required  to  deter- 
mine the  transition  point;  the  transition  may  occur  at  a 
flow  and  RP?!  acceptable  for  the  motor  application.  However, 
the  utilization  of  external,  f orced-convective  cooling  has 
some    major    difficulties   associated    with    it: 

1)       The      source    of     the    external  cooling      fluid    must     be 

easily    accessible,    dependable    (in   a  casualty-control   sense), 

safe,  volumetrically  small,  light,  and  should  be  reasonably 
inexpensive . 
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2)  The  fluid  itself  should  be  of  high  specific  heat, 
light,  inexpensive,  and  sore  importantly,  non-toxic  and 
r.on-f  la  it,  enable. 

3)  The  method  of  transporting  the  fluid  from  the  source, 
through  an  external  heat  exchanger,  through  the  device,  and 
tack,  should  be  extremely  reliable  and  quire  resistant  to 
external  damage.  This  is  an  area  of  serious  consideration 
wher: e  the  reliablility  of  conventional,  rotating  seals  is 
SliS^^^2ilAk2*   -H   lh.2.    least . 

The  alternative  to  the  use  of  external,  forced- 
convective  cooling  is  the  utilization  of  a  closed  device: 
(a)  the  closed,  two-phase  t  her  mosy  phon,  (b)  the  closed, 
rota ting-locp  thermosyphon,  (c)       the   closed   rot ating-loop, 

two-phase  t her mosy phcn,  or  (d)  the  heat  pipe.  These  alter- 
natives remove  the  necessity  of  using  radial  rotating  seals 
in   the    design   of   the    cooling   system. 

The  closed,  two-phase  thermosyphon,  as  shown  in  Figure 
3.6,  must  have  a  sufficient  forced-air  circulation  through 
the  end-tell  regions  to  remove  the  heat.  This  will  reguire 
optimizing  the  length  of  the  condenser  sections,  the  geom- 
etry of  the  fins,  the  number  of  fins  and  the  direction  of 
air  flow  within  the  device  to  maximize  the  heat  transfer 
while  minimizing  the  windage  losses.  The  closed  two-phase 
thermosyphon  is  feasible  at  the  higher  RPM's  that  the  gener- 
ator operates  at;  additionally,  this  method  could  also  be 
feasible  in  the  motcr  application  if  a  heat  pipe  (i.e.,  a 
thermosyphon  whose  walls  are  lined  with  a  capillary  wicking 
material)  is  used  in  lieu  of  the  two-phase  thermosyphon  for 
the    lower   RFH" s. 

The  closed,  rotating-loop  thermosyphon  requires  an  addi- 
tional heat  *  exchange  through  the  shaft  of  the  device, 
possibly  supplemented  with  forced-air  convection  through  the 
hollow   rotor      region    and   by      the    exposed   radial      sections    of 
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the  loop.  This  shaft  heat  exchanger  may  be  in  the  fen  of  a 
rotating  heat-pipe,  as  shown  in  Figures  4.1  and  4.2,  where 
the  external  exchange  is  to  either  ambient  forced-air,  cr  to 
an  external,  closed  water  spray.  The  shaft  heat  exchanger 
may  also  he  cooled  by  forced  convection  using  liquid 
cooling,  with  axial,  packing-type  seals,  external  to  the 
device  as  shown  in  Figure  4.3.  This  closed-loop  thermcsy- 
phon    could    use    either  single-phase   or   two-phase   operation. 

The  closed  two-phase  th ermosyphon  could  also  be  extended 
along  the  shaft  through  a  bearing  plate  arrangement  and 
cooling  air  could  re  supplied  external  to  the  casing, 
simplifying  the  internal  configuration,  as  shown  in  Figure 
4.4. 

The  ability  of  these  devices  to  cool  effectively  has 
been  proven  and  their  usage  in  this  application  is  theoreti- 
cally feasible.  The  exact  correlations  for  the  calculation 
of  the  performance  values  for  these  devices,  especially  in  a 
rotating   reference,    will   require    further  research. 

The  analysis  herein  has  shewn  the  ability  of  th*  cooling 
schemes  discussed  to  remove  the  .requisite  heat  load  in  order 
to  extend  the  life  cf  the  insulation,  increase  the  effi- 
ciency, decrease  the  weight  and  size  of  motors  and  genera- 
tors. This  study  graphically  presents  the  applicable 
correlations  for  the  determination  of  the  heat  transfer 
within  a  rotating,  electric  device.  It  also  discusses  the 
limitations  of  these  correlations.  It  has  been  shown  that 
closed,  two-phase  systems  of  an  inherently  higher  reli- 
ability  are    feasible,    and   deserve   further   evaluation. 
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7*     CONCIDSIONS    ANP   BJCOH MEN DAT IONS 

A.       CCNCIOSIONS 

Eased  upon  the  information  reviewed  and  analyzed  in  this 
thesis,    the   following  conclusions   are   raad^: 

1)  For  ced-convective  liquid  cooling  is  a  practical 
means  for  improving  the  efficiency,  reducing  the  size  and 
weight,  and  extending  the  life  of  high-power  motors  and 
generators.  The  lain  disadvantage  of  forced-convective 
liquid   cooling    is  the   requirement   of    radial  rotating   seals. 

2)  Closed-loop  liquid  thermosy phons  are  feasible  at 
high  EFM.  They  offer  improved  reliability  over  forced- 
convective  liquid  ceding  which  utilizes  rotating  seals. 
However,  this  method  requires  the  use  the  shaft  as  a 
secondary   heat    exchanger. 

3)  The  closed,  two-phase  thermosyphen  is  feasible  at 
the  higher  RPMs  of  the  generators.  This  method  requires 
forced-air  convection  through  the  end-bell  regions  or  the 
extension  of  the  devices  alcng  the  shaft  for  external 
forced-air  convection,  both  of  which  will  add  to  the  windage 
losses. 

U)  For  low  RPM  applications,  such  as  motors,  the  use  of 
a  heat  pipe  will  be  required  to  overcome  the  inability  of 
the  lew  acceleration  field  to  transport  the  cooling  fluid. 
Both  the  closed,  two-phase  thermosyphon  and  the  heat-pipe 
require  forced  convection  at  the  condenser  ends  to  remove 
the  heat  either  internal,  or  external  to  the  device.  This 
may    add   to    the    windace   losses. 
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5)  Correlations  and  theoretical  analyses  are  net  avail- 
able   in    the   literature   fcr: 

a)  the    counter-current  entrainment    limit    for   two- 
phase    t her mesyphens    under  rotation   and 

b)  the      transition    from   laminar   tc    turbulent    flow 
for    rotating  systems. 

E.       RECOMMENDATIONS 

The  following  recommendations  to  continue  this  work  are 
made: 

1)  Testing  of  the  f crced-convective  liquid  cooling 
scheme  should  be  accomplished  using  the  DTNSHDC  test  rig  to 
confirm  the  analysis  and  the  validity  of  the  available 
correlations,  or  to  obtain  proper  correlations  for  the 
specific   geometry  and   utilization   being   projected. 

2)  An  experimental  program  should  be  devised  and 
conducted  tc  determine  the  entrainment  limitation  in  a  hori- 
zontal, heated,  channel  rotated  about  a  parallel  axis  with 
counter-current    flow    cf   liquid   and   vapor. 

3)  An  experimental  program  should  be  developed  to 
determine  the  optimum  condenser  length  and  proper  fin 
arrangment  for  the  use  of  both  heat  pipe  and  closed  two- 
phase  ther mosyphon  systems  for  rotor  cooling.  Consideration 
cf  the  structural  integrity  must  be  included  for  the 
condenser   end-lengths. 

4)  An  experimental  program  should  be  developed  for 
correlating  the  laminar- to-turbulent  transition  during 
f orced-convect ive  liquid  cooling  in  a  heated,  horizontal, 
rotating  pipe. 

5)  An  experimental  program  should  be  developed  to 
modify   the      inside    gecmetry    of   the      closed-loop   thermosyphon 
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(and    axial,    closed,    two-phase   thermos yphon)      to   optimiz- 

heat-transfer  and  the  friction    losses    within   the    device 


li  e 


6)  A  detailed  analysis  should  be  made  of  the  shaft 
cooling  potential  fcr  the  rotating-loop  thermosyphon  with 
regard  to  the  heat-transfer  enhancement  cf  the  shaft  to  the 
fluid    and  the   mechanical   design   of   the    shaft. 

7)  An  experimental  program  should  be  developed  to 
analyze  *:he  closed,  two-phase  thermosyphon  with  radial 
condenser  sections  and  the  closed- loop  two-phase  thermosy- 
phon. These  two  devices  may  experience  start-up  balance 
problems. 

8)  Various  configurations  should  be  analyzed  for  their 
potential  for  heat  transfer  to  include  a  combination  of  the 
systems  described  herein,  as  depicted  in  Figures  4.1  through 
U.4. 

9)  The  concept  cf  modifying  the  geometry  of  the  rotor 
windings  within  these  devices  in  order  to  minimize  heat 
conduction  path  lengths  to  the  casing,  shaft,  etc.,  and 
thereby  the  internal  temperature  gradient,  should  also  be 
investigated. 
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APPODIX    A 
GOVERNING    EQUATIONS    FOR   LAMINAR    CONVECTION    IN    UNIFORMLY 
HEATED   HORIZONTAL    PIPES    AT    LOW    RAYLEIGH    NUMBERS 

The      following    development,        by      Morton   [Ref.    11]*.        is 
presented      herein      fcr   completeness.  The      steady      laminar 

motion  of  a  fluid  in  a  horizontal  circular  pipe  of  radius  a, 
the  walls  of  which  are  heated  uniformly  so  that  a  constant 
temperature  gradient  t  is  maintained  in  the  direction  of 
the  axis.  The  flow  will  be  referred  to  in  cylindrical  coor- 
dinates (r,^i  rz)  with  <j>  measured  from  the  upward  vertical 
and  the  z-axis  along  the  axis  of  the  pipe;  the  velocity 
components  are  denoted  by  u,v,  and  w.  The  effects  of  dissi- 
pation and  cf  the  pressure  term  in  the  energy  eguation  will 
fce  neglected  and  that  variations  in  the  density  due  to  the 
temperature   differences   are      so    small    that    they      only   effect 

the      buoyancy      term.        Thermal      conductivity      and      kinematic 

i 
viscosity   are   assumed     to   be   constant    (which      will    introduce 

guantitative  errors    into  the   solution,    but    should    not   change 

its    general   character). 

The   eguation   of    continuity    is   then: 

3ru/3r    ♦    3v/3  4>    ♦   3rw/3z    =0,  (A.1) 

and   the   energy    eguation   is: 

u(3T/3r)     ♦   v/r(3T/3$)    +    w  (3  T/3  z)     =    k   72T,  (A. 2) 

where  V2=32/3r2  ♦  VrO/Sr)  +  1/r2  (3  2/3cj>  2)  +  32/3z2  and  T  is 
the  local  fluid  temperature.  The  momentum  equations  can  be 
written    in    the    form    cf   equations    (A. 3)     through    (A. 5). 
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a(3u/3r)  *Vr(3u/3<j>)+w(3u/3z) -v2/r   =   -1/p(3p/3r) 

v(    72u-u/r2-2/r2  iSv/d-i)-    a  (T   -T)  cose    ,  (A.  3) 

w 

u(3v/3r)  +  v/r  (  3v/3<J>)  +w  (  3v/ 3z)  =uv/r   =   -1/pr  (3p/3<j)) 

v{    72v*2/r2  (3u/3<j>)-v/r2)+    g(T  -T)  sin<|>    ,  (A. 4) 

W 

u(3w/3r)  *v/r  (  3v/3o) +w(3w/3z)     = 

-1/P(3p/3z)+v    v2w.  (A. 5) 

In   equations    (A.  3)       and    (A.  4)  ,         the    buoyancy    force    has   been 

calculated    relative    to   the    fluid      at   the   same   level   adjacent 

to    the   pipe   wall   and    ^he  remaining   distribution   of    force   has 

been    absorbed   into    the   pressure,    p. 

For     steady    convection      sufficiently   far      from   the      pipe 

opening     to     avoid    inlet -length      effects,        the      temperature 

throughou-1-    the      pipe    increases      uniformly    with     the    distance 

along    the    pipe    axis.        Hence,        the   distribution    of    buoyancy 

force    in   sections      of  the  pipe    is   independent   of      z;      as   the 

secondary   flow    is      caused   by  the    buoyancy      forces,      the   flow 

field   must    also    be    independent   of   z.        It    follows    that   there 

should    be      a   similarity      solution    with      u,v,w,      and      T    -T   as 

w 

functions   of  r    and    $      only,       and    with    p=Yz  +  P(r,  <t>   )        (where   P 
contains   the  terms    absorbed    into    p)  . 

The   continuity    equation    reduces    to   the    form: 

3(ru)  /3r   ♦    3v/3<J>  =    0; 

hence,      a      dimensionless  Stokes      stream    function    ty        can     be 
introduced    in   such    a    way  that: 

ru/v=3^/3cf),  v/v  =  -  3i|j/3r. 
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The    main  equations    (A.1)    through    (A.  5)     may   be    reduced    to 

non-dimensional    form    by   the    transformations   r=aR,      w=(v/a)W, 

and    T   -T=iaPr    ,      where   Pr  is      the   Prandtl   number     (v/a)  .        If 
w 

the      pressure   is      eliminated      between      equations     (A. 3)         and 
(A. 4),    the   resulting    momentum   equations   are: 

Pa{(38/3R)  sin(j>      ♦     1/R  (3 9/3 cj>)  cos   <£}  r  (A. 6) 

7*W*    1/BC(3^/3R)  O/ad))- 

(3i|;/3<t>)  (3/3R)}  K    +  4Re    =0  r  (A. 7) 

and    where, 

Vf   =         32/3R2    +    1/R(3/3R)  ♦    1/R2(    32/9^     z)  . 

The  Bayleigh  number,  Pa,  is  given  by  gSia*/(av  )  and  the 
Reynolds  number  is  the  normal  Reynolds  number.  Re,  for 
Foiseuille  pipe  flow  based  on  the  pipe  diameter  and  the  mean 
velocity  across    a   pipe  section    is    (a/v)  (  Ya2/4U    )• 

The  solution  of  the  energy  equation  depends  on  the  temp- 
erature en  the  boundary  condition  at  the  wall  of  the  pipe. 
Pipes  will  usually  have  reasonably-thick  walls  of  material 
with  thermal  conductivity  much  higher  than  that  of  the 
fluid,  so  there  will  be  little  variation  in  wall  tempera- 
ture around  the  pipe  circumference.  This  will  be  specially 
so  for  slow  rates  of  heating  when  asymmetries  of  the  flow 
will  have  small  amplitude.  Hence,  write  T  =T0+  z,  where  TO 
is  the  wall  temperature  in  the  section  containing  the 
origin.  Although  there  is  uniform  heat  transfer  per  unit 
length  of  pipe,  the  local  heat  transfer  will  be  slightly 
greater  near  the  bottom  of  the  pipe  than  near  the  top.  With 
this  assumption,  equation  (A. 2)  rsduces  to  the  non- 
dimensional   form   of    equation    (A. 8). 
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72e    +  (cr/fij  [(3ip/3B)  (3  6/3*)  -(3.1-/3*)  (36/3r)}+w=0  (A.  8) 

The   boundary  conditions    are: 

u,v,w,    and    3      are    zero    on    r=a; 
u,v,w,   and    9      are    finite    on    r=0. 

(A. 9) 

Although  solution  of  equations  (A. 6)  through  (A. 8),  satis- 
fying conditions  of  (A. 9)  is  considerably  difficult,  succes- 
sive approximations  to  the  solution  can  be  obtained  by 
expanding  ip,  w,  and  9  as  power  series  in  the  Rayleigh 
number,  provided  that  this  is  numerically  small.  Supposing 
that: 

^=Ra^  ^Ra2^   +. . ., 

W=W  0+SaW  ^Ra  2W2  +  .  .  .  , 

9=90  +  pa9  i+Ra262  +  .- •    . 

(A.  10) 

The    leading     term    ^ 0    of      ty      must    vanish    because   there      is   no 

circulaticn    when      A    is   zero,      but    90      is   not   zero      since   the 

difference    in      temperature    between   a      fluid  element      and   the 

neighboring   wall      is    proportional      to    t9    •  When    relations 

(A.  10)    are    substituted   in   equations    (A.  6),     (A. 7)  ,    and     (A.  8)  , 

three   sets   of   equations    for    the      funtions    \b     ,      w     ,      and    9. 

i  1  1 

are    obtained  by    equating  coefficients    of    powers   of   Ra. 

From  equation    (A. 7),    the   basic   equation   is: 
72ro    ♦    URe  =    0,  (A.  11) 

which    has   a    solution: 


W  Q=    Re(1-r2)  (A.  12) 

which    satisfies      the    conditions   W  0=0      at   R=1,      »0      finite   at 
R=0.  This  is      ordinary   Poiseulle      flow      under   a      pressure 

gradient   -Up^Re/a3    in  an   unheated   pipe. 
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The    first   of     the   equations   derived    from      equation     (A. 8) 
for    the   temperature    distribution    is: 


-2e0+*  =  o, 


(A.  13) 


subject    to    90=0    at    3=1,  80    finite   at   R=0,      and   the    solution 

is      axialiy   symmetrical  as      no   account      has      been   taken     of 

gravity    at    this    stage  of  the   approximation.        This   is   satis- 
fied   by: 


60    =(V16)  Ee(1-R2)   (3-R2)  , 


(A. 14) 


which    is   the  customary   solution    for    forced    convection. 

The    first-order    approximation   for    \l>     satisfies    the   equa- 
tion: 


^^=090/33)  sin    $ 


(A.  15) 


obtained  as  the  coefficient  of  Ra  from  equation  A.  6,  and  the 
boundary  conditions  a^j/SR,  9^1/3$=0  on  R=1  ,  and  R-M9^i/3^)» 
34>i/3  R  remain  finite  at  R=0.  If  ty-  is  assumed  to  have  the 
form  <p  ,  (R)  sin<+>  ,  the  dependence  of  equation  A.  15  on  <b  is 
eliminated,    and    ip    (R)    is  easily   found: 


ty    =-  (1/U608)  ReR(  1-R2)  2(1  0-R2)  sin  (J) 


(A.  16) 


Similarly,      the    first-order    approximation    to    w^satisfies   the 
equation : 


72W    =  (1/R)  (3^x/3<J))   (3W./3R)  , 
l     l  u 


(A.  17) 


and  the  boundary  conditions,  Wx=0  at  R=1  and  wx  finite  at 
R=0.  The  dependence  of  equation  A.  17  on  <£>  is  eliminated  by 
talcing    W1=W1  (R)  cos^    ,   whence 


102 


W  j  =  -(  1/18U32  0)  Re*R  (1-R2)  (49-5  1R2+ 19R*-R&)  ccsi  (A.  18) 

The  first-order  approximar  icn  to  0  satisfies  the  equa- 
tion: 

V29    =(o/R)  (3<K/3<J>)  <39n/3R>  *H,  '  <A«  19) 

and  the  boundary  conditions,  9  =0  at  R=1  ar.d  6  finite  at 
B=0.      Hence, 

9  !=-  (1/22118400)  Re2S(l-R2)  { ( 38 1 + 1 325Q) - ( 354 +16 7  5a ) R* + 
+  (146+925o    )  R*-  (  29+200a)  R  6  +  (  1-2a)  R8}  ccs^     .  (A. 20) 

This  completes  the  approximate  solution  to  the  first 
order  in  Ra  for  convection  in  a  horizontal  pipe,  which  is 
heated  so  that  there  is  a  constant  temperature  gradient 
along  the  walls  and  uniform  temperature  around  the  girth. 
It    may      be    noted      f rem  equations      (A.  10)  ,       (A.  12)  ,  (A. 14), 

(A.  16),  (A.  18),  and  (A.  20)  that  the  full  convection  solution 
depends  essentially  en  the  product  RaRe.  This  is  clear  if 
it  is  recalled  that  Ra  is  proportional  to  the  increase  in 
wall  temperature  along  the  pipe  and  R2  is  proportional  to  a 
characteristic    velocity     of    the      flow   along      the    pipe.  An 

increase  in  Ra  means  that  the  fluid  will  be  carried  through 
a  larger  temperature  differences,  and,  hence,  there  will  be 
increased  buoyancy  forces;  but  the  same  effect  can  be 
produced  by  an  increase  in  Re,  and  hence  in  the  velocity  of 
the    main    flow. 

The  two  most  interesting  consequences  of  such  a  solution 
are  the  modification  cf  the  heat  transfer  due  to  secondary 
flow,  and  the  effect  of  heating  on  the  tangential  stress  at 
the  wall  (or  on  the  rate  of  flow  for  a  given  pressure 
gradient).         However,      as      both    W  x  and     9  vary    with   cos  cf  » 
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neither  the  flux  ncr  the  heat  transfer  across  the  whole 
section  of  the  pipe  will  be  changed  by  the  first-cr3er 
approximation   although  there    will    be    local    variations. 

To  find  *-he  changes  in  total  flux  and  heat  transfer,  it 
is  necessary  to  proceed  to  the  second-order  approximations, 
and   the   revelant   equations    for   these   are: 

~*    i,  =(VR)   (3<     v2>i;-  ,xb  -)/3(R,<!>    }} 

+  (36    /3    R)  Sin *♦  (1/R)  ( 9  6/3 cj))  cos*    , 
i  r  i 

V2W    =(1/P.)   (3(1    ,*    )/3  (R,  <D)}+(VR)  0(Wn.,^  J/3(R,<J>)}  , 

12  l  *  U  ^ 


=  (a/R)  {3(9    ,110/9  (R,4>)  }*(a/R)  OO^'Jj/a  (R,4>)  } -*2  - 

(A. 21) 


1      2 


The  only  difficulty  in  solving  these  equations  is  numerical 
tediom  and  the  remaining  second-order  solutions  are 
presented   in  [Ref.    11]. 
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&PPJJDIX    3 

SINGLE-FHASE    ANALYSIS    (COMPUTER    PROGRAM    AND    RESULTS) 


The  following  program  was  written  to  give  the  results  of 
the  various  correlations  for  the  heat  transfer  for  th« 
forced,  single-phase  convection  in  the  rotating  reference  as 
disscussed  in  Chapter  II.  The  program  is  so  written  that 
modification  to  include  other  correlations  is  guite  easy. 
It  was  written  in  HP-Basic  for  the  HP-9826  computer. 
Sub-programs  for  the      graphics  are   not    shown.  Sample   runs 

are    included  immediately  following   the   listing. 

-Pages     105-109      Program    Listing 

-Pages    110-111       Thermophysical    Properties    (Common   to 

Appendices   B    through   D)  . 

-Page       112  Sample    Calculation 
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!  rI_E     lAME:  -,  ■■ 

'  CREATED:  .annarv    13.    _ 934 

lu!  urtS  '    REVlSrfj;    -  ^ r' >■  i^.a r v    -  ~.       >  -.-.« 

U'JU! 

040  COM  /Cmin/  Xmin. Ymin.Sf x .Sfy 

050  PRir-TE-  IS  1 

360  BEEF 

0^0  PRIN1  JSING  "2X. ""Defau.lt  values:  

030  D=4.763E~3     !  Tune  o.ameter    (m) 

'090  L=.8230        !  Tube  length      (m) 

1100  R=.402        !  Radius  of  rotor  <n> 

niO  PRINT  JSING  "4X.""Tube  diameter   =  "" .Z.4DE ,""  (m) :D 

1120  PRINT  USING  "4X  .  ""Tune  length     =  M,,.DD.D.""  <m) ;L 

I  130  PRINT  JSING  "4X. ""Radius  of  rotor  =  "".Z.3D,""  <m) :R 

1  M0  BEEP 

1'50  INPUT  "OK  TO  ACCEPT  DEFAULT  VALUES  < 1 =Y . 0-N>?".Id 

1160  BEE=> 

••70  INPUT  "LIKE  A  HARD  COPY  ( ! -Y .0=N>?" . Ihc 

1180  IF  Ihc=l  THEN 

1 '90  PRINTER  IS  701 

120C  E_SE 

1210  -RINTER  TS  ! 

1220  END   I?' 

1230  Ir    Id-1    THEN    1 300 

1240  BEEP 

^250  INPUT    "ENTER    "TUBE    DIAMETER    <m)".D 

1 °S0  BEEP 

•270  INPUT  "ENTER  TUBE  LENGTH  (m)",L 

1280  BEED 

'290  INPUT  "ENTER  RADIUS  OF  ROTOR  (m)".R 

'300  PRINT  USING  "10X.""*»*  Geometric  Variables  »**  

'310  PRINT  :SING  "14X.""Tube  diameter  <D>  --    ""  .Z.3DE .""  <n) :D 

"320  PRINT  USING  "14X.""Tube  ienqtn  (L)    -  "".DD.DD.""  <m> :L 

:  330  DRINT  USING  "14X.""L/D  =  ""  .5D.D":!_/D 

'340  PRINT  USING  "MX, ""Rotor  radius  (R)   =  "".Z.3D.""  Cm) ;R 

■350  PRINT 

•360  °RINT  'JSING  "10X.""ComDutec  ,a:ues  are:  

1370  PRINTER  IS  1 

•380  BEEP 

'390  DRIN~  .SING  "2X. ""SELECT  OPTION: 

MOO  PRINT  USING  "4X.""1   Sinale  point 

1410  PRINT  JSING  "4X.""2   NU  versus  RPM 

1420  PRINT  USING  "4X.""3   Nu  versus  Re 

1430  DRINT  USING  "4X,""4   Nu  versus  Prod 

1440  INPUT  Op 

1450  IF  0p>3  THEN 

:460  PRINT  USING  "4X.  ""SELECT  OPTION: 

1470  PRINT  USING  "4X.,,,,2   VARY  RPK 

1480  PRINT  USING  "4X,""3   VARY  Re 

1490  INPUT  Opx 

'500  IF  0dx=2  "«en  up=2 

1510  IF  0px=3  THEN  0p=3 

1520  END  IF 

1530  IF    0p>1    THEN 

1540  BEEP 

1550  Jj=0 

1 5B0  PRINTER    IS    1 

'570  PRINT    USING    "4X . ""SELECT    OPTION: 

1580  PRINT    USING    "4X.""      1  =Di  ttus-Boel  ter  .    2=Nakayama.    3  =  Nakayama   b    Fuzioka' 
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1t>90  PRINT  USING  "4X,"'   4=Woods  ~t   Morns.  5=Stephenson.  6=Woods&Morns*'2 ,  7=Lam 

i  nar  .  o  =  f-<a>.  a.v  ina-  1 " ' 

1 600  INPUT  Ian 

I62G  mruT  'i_I\E  Tu  ri_ui  i  i"  i  .  O^N)'"  .Ukpiot 

1630  IF  Oirplot-1  THEN 

1640  CALL  Piot(Jj.oC.jd.X.Y.Type.Cx) 

1650  BEEP 

'560  INPUT  "SOLID-1 ,DASH-DflSH-2.DQT-D0T-3".Tvpe 

1670  END  IF 

'580  END  IF 

1690  Nstep-100 

'"00  IF  Dp=2  THEN 

1710  BEEP 

'720  INPUT  "ENTER  RPM  RANGE  < MIN  .MAX )" .Rpml  .Rpmh 

1730  Rpm=Rpml 

1740  ELSE 

1750  BEEP 

1760  INPUT  "ENTER  RPM  OF  ROTOR". Rpm 

1770  END  IF 

! 780  BEEP 

1790  INPUT  "ENTER  OPTION  (1=Q,2=Mf ,3=Dt>"  .  Io 

1800  IF  Ihc-1  THEN  PRINTER  IS  701 

1310  IF  lhc=0  THEN  PRINTER  IS  1 

1820  PRINT 

1330  DRINT  USING  "10X,""**«  Operating  Variables  ***  

1840  PRINT 

1350  PRINT  USING  "1  OX . ""Input  variables  are:  

I860.  IF  Io=1  THEN 

!870  3EEP 

1380  INPUT  "ENTER  COOLANT  MASS  PLOW  RATE  (kg/s>".Mf 

1890  BEE* 

'900  °RINT  USING  "MX. ""Coolant  nass  Flow  rate  =  "".Z.3DE,""  (kg/s> ;Mf 

1910  INPUT  "ENTER  INLET  AND  OUTLET  TEMPS  <C>".Ti.To 

1320  3RINT  USING  " MX . ""Coolant  inlet  temp     -  "".DD.DD.""  <C) :Ti 

1930  PRINT  USING  "1 4X . ""Coolant  outlet  tenc    =  "".DD.DD.""  (C) :To 

1940  END  IF 

'950  Ir  Io=2  THEN 

1960  3EEP 

1970  INPUT  "ENTER  HEAT  LOAD  <W)".Q 

1980  PRINT  USING  "14X.""heat  load  =  "".4D.DD.""  (W) ;Q 

'990  BEEP 

2000  INPUT  "ENTER  INLET  AND  OUTLET  TEMPS  (C>".Ti.To 

2010  PRINT  USING  "1 4X . ""Coolant  inlet  temo     =  "".DD.DD.""  <C> :Ti 

2020  PRINT  USING  "14X . ""Coolant  outlet  temp     ■  "".DD.DD,""  (C> :To 

2030  END  IF 

2040  IF  Io»3  THEN 

2050  IF  Ut>0  THEN  2200 

2060  BEEP" 

2070  INPUT  "ENTER  HEAT  LOAD  <U)".0 

2080  PRINT  USING  "14X.""Heat  load  =  "".4D.DD,""  (W) :0 

2090  BEEP 

2100  IF  0p<3  THEN 

2110  INPUT  "ENTER  COOLANT  MASS  rL0W  RATE  <kg/s)",Mf 

2120  PRINT  USING  "1 4X . ""Coolant  mass  flou  rate  =  "".Z.4DE.""  <kg/s> ;Mf 

2130  ELSE 

2140  INPUT  "ENTER  MASS  FLOW  RATES  (MIN .MAX)" ,Mf 1 .Mf h 

2150  Mf-Mfl 

2'S0  END  IF 

2170  3EED 

2^80  INPUT  "ENTER  COOLANT  INLET  TEMP  (CV.Ti 
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:-3j 

22;0 
222u 

2230 

2240 
2250 
22S0 
2270 
2230 

2300 

2310 

2320 

2230 

2340 

2250 

2360 

2370 

2380 

2390 

2400 

2*10 

2420 

2430 

"  :  ~o 

2440 

2450 

2460 

2470 

7a 

2^80 

2490 

2500 

2510 

2520 

2530 

2540 

2550 

2560 

2570 

2580 

2590 

2500 

2510 

2520! 

2530 ! 

2640 

2550 

2660 

2570 

2580' 

2590 

2700 

2710 

2720 

2730 

2740 

2750 

2760 

2770 


PRINT  USING  "•4X.,,,,Cooiar.t  inlet  temp     =  ",',DL.DC, 

■    if*~_ 

id.-v  iot7i)*,5 

Lp=hNUpu<  Ta> 

Mu-FNMuu<Ta) 

<  =  FNKL«(Ta) 

Pr«FNPru(Ta) 

Rho=FNRhou<Ta> 

Nuv=Mu/Rho 

3=Mf«CD*(To-Ti) 

PRINT  USING  ,,14X,,,MHeat  load  <Q>  =  ,,".4D.D.' 

END  IF 

IF  Io=2  THEN 

Mf =Q/(Cd*(T0-ti ) > 

PRINT  USING  "14X. ""Coolant  mass  flow  rate  =  "".Z.4DE, 

END  IF 

IF  Io=3  THEN 

Toc  =  Ti+Q/(!^f*CD) 

IF  PBS(To-Toc>>.01  THEN 

To=( To+Toc>* .5 

GOTO  2220 

END  IF 

IF  Jj=0  THEN  PRINT  USING  "1 4X . ""Coolant  outlet  temp 


( C ) ' 


<W)' 


(kg/s) 


;Mf 


=  "".3D.DD,""  (C)' 


END  IF 

IF  Jt=0  THEN 

PRINT 

PRINT  USING  "10X. ""Fluid  properties  evaluated  at  "".DD.DD.  ""    (C)  are:"" 


DRINT  USING  "14X. ""Specif ic  heat  (Cp) 
PRINT  USING  "MX."" Viscosity  (Mu) 

PRINT  USING  "MX  .""T-.ermal  cone!  <k> 

DRINT  USING  "14X.""Prandti  numoer  (Pr) 

PRINT  USING  "14X.  ""Density  <Rho> 

PRINT  USING  "14X. ""Kinematic  vis  (Nuv) 

3eta=[rNBeta(Ta) 

PRINT  USING  "14X.""Coef  ther  exp  (Beta)    =  " 

PRINT 

PRINT  USING  "10X.""*»*  Calculations  *** 

PRINT 

PRINT  USING  "10X.""'Drel  iminary  calculations:' 

END  IF 

Re=4*r1f/(PI»D*Mu) 

Friction  Factors  for  stationary  reference 

IF  Re<2*1 .E+4  THEN  2670 

r=. !34/Re" .2   !I4D  Eqn  8.21 

GOTO  2590 

F=.316/Re".25   ! IdD  Eqn  8.20 

Qpd=Q/(PI*D*L> 
Qp=Qpp*PI-D 
Vm«4*Hf/<Rho»PI»D~2> 
Dps-F»Rho»V«T2*L/<2»D) 

0meaa=Rpm*2»PI/60 

Ro=0m/(0mega*D) 

Jay=Omega*D '2/Nuv       !"rom    STEPH-Rotat lonal    Re 

KiL'  =  FNKu(Ta> 

Dtt-<To-Ti)/L 


4D.D.""    <J/kg.K> :Cp 

Z.4DE.""    (N.s/m'2) :Mu 

Z.4D.""    (W/m.K) :K 

,Z.3D":Pr 
4D.D.""    (kg/m    3) :Rho 

,2  .4DE  ,""    (m    2/s> :Nuv 

Z.4DE 


"    (1/K) :Beta 
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"SO 
2790 
2800 

*>  •".  *     ~s 

CO   t  \J 

282U! 

2830! 

2840 

2850! 

2860! 

2370 

2320! 

2890! 

2900 

er 

2910! 

~1820 ! 

2930 ' 

blent 

2940 

2950! 

2960! 

2970 

2980 

2990 

3000 

3010! 

2020! 

3030 

2040 

3050! 

2060! 

3070 

3080 

2n90' 

3100! 

3110 

3!20 

3130! 

~:40! 

3150 

3'60 

3!70 

3'80 

3<90 

3200 

2210 

3220 

3230 

3240 

2250 

2260 

3270 

3280 

3290 

2300 

3310 

3320 

3330 

3?40 

3250 

3360 


Gr=R* ( Omega " 2 ) *Be ta*  <  ( D/2 )  4  >  * L!  t :  / ' I uv  2 

Ra=Gr»Pr 

Ganma=(Re ' (22/1 3>>*((Gr*(Pr  .6> >"<- 12/13)) 

T     ■"»    /  /  f ■    ~\      C   \ 

Laminar  Value 
Nulam=48/1  1 

Turbulent  Value  from  Di  t  tus-Boe  1  ter 
Nudb=.023*Re~.8*Pr'  .4 

Steonenson   Correlation 

Nus=<Pr    .4/ .367>*.007>Re~ .88*Jay    .023    !fiir ,Turb-EQN2.23-corrected    for   uat 


Nakayama  Correlation    for    turbulent    conditions 

Nuna*c  =  (Re~  .8)»<Pr"  .4>*(Tg'  (1/30)  )*<  1  +  .  01  4/(  Tg)    <1/6>)«.033       !    Eqn'< 

Ty  =  Nunak/(Re" .3*Pr~ .4) 

Naxayama   Correlation    for    laminar    conditions 

J2=J*2 

Sq-2/1  W!+SQR<  l+77/4*<  |/Pr"2>>) 

Cf=1-.486*((3*Sq-1)' .4/(Sq*(Sq*Pr*SQR(5)+2))>»J2/(Ra«Re)     .6 

Nunak 1=48/11*. 1 91 /Sa*(3*Sq- 1 ) " .2*(Ra-Re)  "  .2*1 /( I+( ! /10*Sq*Pr) > 

Nakayama/Fuz  loka  Correlation    for    Radial    Pines 
Nunf-(.014/.023>-(Re/Ro'2.5>' .124»Nudb      !TURB-EQN   2.21 
Mus  =  Mu    !  Temp   Value   Mus'S>Mu 

Sieaer-Tate  Correlation  for  Turbulent 
Nust=.027*Re" .8*Pr' .3333»<Mu/Hus>'  . 14 
Prod=Ra*Re*Pr 

rtoods-Morris  Correlation  for  Laminar 
Nuum=.262*Prod' .173*48/1  1.  !  Eqn  2.22-WM3 

J:=jay/8 


>.9  (Tur 


rtooas-Morris  Correlation  for  Racial 
Nuwm2=.C15*Re*  .78»JT  .25  !from  um2 
IF  Jj-0  THEN 


ipes 


""Reynolds  nujnoer  (Re) 
""Friction  factor  <stat> 
""Heat  flux  <Qpp> 
""Mean  fluid  ve 1  (  Vm ) 

""RPM 


••".Z.3DE" 

:Re 

"".Z.3DE" 

:  r 

"".Z.3DE. 

""    (W/m'2) 

"".Z.2DE. 

""    (m/s) :V 

"".4D.DD, 

"  :  Rpm 

DRINT  USING  "MX 

PRINT  USING  "14X 

°RINT  USING  "14X 

PRINT  USING  "14X 

PRINT  USING  "14X 

END  IF 

IF  lhc=0  THEN  PRINTER  IS  1 

PRINT 

IF  Jj-0  THEN  PRINT  USING  " 1  OX . ""Resul ts :  "" 

IF   Okplot=0    OR    <0kplot=1    AND   NstepOl)    THEN 

PRINT   USING    "10X.""  Nudb  Nunak         Nunf 

PRINT    USING   "13X.6(3D.DD.2X)";Nudb. Nunak .Nunf .Nuwm .Nus .Nuwm2 

PRINT   USING   "14X.""        Ra  Ro  Ra*Re*Pr  Mf 

PRINT   USING   "14X.4(Z.2DE.2X) . 3D .DD ,DD .DDD . Z . DD" :Ra.Ro .Prod ,Mf . 

PRINT   USING   "MX.""        RaRe  Tg  Ty 

PRINT   USING   "13X.D.2DE.2X.Z.^DE.2X.Z.DDD":Ra*Re.Tg.Ty 

END    IF 

IF   Okplot-1    THEN 

IF    Ian=1    THEN    Y=Nudb 

IF    Ian=2    THEN    Y-Nunak 


;Qpp 


Nuwm 


Nus 


Nuuim2' 

Mak  1    "• 
Nunak  1 
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3270 

IF  ian=3  THEN  Y=Nu.nf 

3380 

IF  Iar>=4  THEN  Y  =  Nuwm 

3390 

IF  Ian=5  THEN  Y  =  Nus 

"'  a  0  n 

7-  T3"  =  6  TlJE^  v='" <*<'■' 

3410 

IF  Ian  =  7  THEN  Y=Nulam 

3420 

IF  Ian=8  THEN  Y  =  Ntir,ak  1 

3430 

IF  Qp=2  AND  0px<>2  THEN 

X  =  Rpm 

3440 

IF  0o=3  AND  0dx<>3  THEN 

X=W 

3450 

IF  Qpx-2  GR  0px=3  THEN  X-Prod 

3460 

J  j  =  J  j  + 1 

3470 

CALL  Plot(Jj.Jc,Jd.X.Y." 

Fype  .  Cx  > 

3480 

END  IF 

3490 

IF  Qp=2  THEN 

3500 

Rpm=Rpm«10~<Cx/NsteD> 

3510 

IF  Rp(n>Rpmh  THEN  3590 

3520 

GOTO  2730 

3530 

END  IF 

3540 

IF  0p=3  THEM 

3550 

Mf-Mf»10~<Cx/Nstep) 

3560 

IF  Mf>Mfh  THEN  3590 

3570 

GOTO  2040 

3580 

END  IF 

3590 

BEEP 

3600 

INPUT  "ANOTHER  RUN  (1=Y 

,0-N)?",Ir 

36'  0 

PRINT  "PU" 

3620 

IF  Ir=l  THEN  1250 

3630 

INPUT  "WANT  TO  LABEL?<1' 

■Y.0-N>",I1 

3640 

IF  11  =  1  THEN  CALL  L.aoei 

3650 

END 
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01  c 

: 


'040 
1050 
1060 
1070 
I  080 
1090 
1  '00 
1110 
'20 
130 
'40 
1 50 
•50 
1170 
1  180 
i190 
:200 
1210 
•220 
1230 
!240 
'250 
•2S0 
1270 
•230 
'290 
•200 
1310 
•320 
1330 
•  340 
'350 
'  380 
1370 
'380 
1390 
'400 
1410 
•420 
•430 
^440 
1450 
1460 
1470 
'480 
•490 
1500 
'510 
1520 
'530 
1540 
1550 
•580 
1570 
1530 


-  /!   ^'•~-oc;tr.  ,"i, 


)/(K <7>*( 1-~) *2+K<8> > 


DEF  FNPvst(~steam) 
olf   K<8> 

'I'd     _7     -  -  '  1  "'-i^D  J  ,  -  i1-      iSiTJ^CQC     -iCfi     nncc;^) 

un i n   s . i o/ 1 1 732 .20.575uG7u. i  .  ES  •  o 

READ   K<~> 

T=(Tstean+273. 15>/847 .3 

Sujti  =  0 

FOR   M=0    TO   4 

Sum=Sum+K(N>*<  ]-J)    (N^l  ) 

NEXT   N 

8r  =  So/n/(  !*<  l+K<5)*<  !-T)+KC6>«U-T>~2>  >-<1~ 

cr=EXP(Br) 

P«22120000*Pr 

RETURN  P 

FNEND 

DEF  FNHfa(T) 

Hfg-2477200-2450»(T-10) 

RETURN  Hfc 

FNEND 

DEF  FNMaui(T) 

A=247.3/(T+133. 15) 

Mu=2.4E-5*10"fi 

RETURN  rtu 

FNEND 

DEF  FNVvst(Tt) 

P=FNPvst<Tt) 

T=Tt+273.15 

X=1500/T 

F1  =  1 /(  t+T«1  .E-4) 

F2  =  (1-EXP(-X))~2.5*EXP(X)/X'  .5 
8= . 00 ! 5*F • - . 000942*F2- . 000«882*X 
K=2»P/(461 .52*T> 
V-<1  +  (1  +2*8*0  '.5)/K 

RETURN  V 

FNEND 

DEr   FNCouifT) 

Cpw-4 .21 i20858-T«(2.26o2SE-3-T*<4.42361E-5+2.7 1423E-7*T> ) 

RETURN  Cdu«*1000 

FNEND 

DEF  FNRhou(T) 

Ro=999. 529*6+  T*<  .  0  12S9-T-<5  .48251  3E-3-T*  1  .234i47E-5>  > 

RETURN  Ro 

FNEND 

DEF  FNPrui(T) 

Pnj  =  FNCpw<  T)»FNMaui<  T  )  /FNKut  T) 

RETURN    Pru 

FNEND 

DEF    FNKui(T) 

X=(T+273.15)/273.15 

<u»=-.9224  7  +  X*(2.8395-X*(1  .8007-X»(  .52577- . 07344»X) >  ) 

RETURN   Ku 

"NEND 

DEF   FNTanh(X) 

P=EXP(X) 

Q=EXP(-X) 

7ann=(P-Q>/<P+Q> 

RETURN   Tanh 

FNEND 

DEF    FNHf(T) 


.cc,- 
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'  590  -if  =  7«<4.  203849-T«<  5.38  T  32E-  4-  T*4.55  1F03'7E-S>  ) 

1600  RETURN  Hf »1 000 

1610  FNEND 

1P20  DEr  FNTw-;^^ 

'  oZ»j  i  a-  i  1  0 

1640  T1=10 

1550  Ta=<Tu+Tl)*.5 

1660  pc=FNPvst<Ta) 

T  570  Ir  ABS<(P-Pc>/P>>.001  THEN 

1S80  IF  Pc<P  THEN  Tl=Ta 

'690  IF  Dc>D  THEN  Tu=Ta 

1700  GOTO  1650 

1710  END  Ir 

1720  RETURN  Ta 

1730  FNEND 

1740  DEF  FNBeta(T) 

1750  Rop=FNRhou(T+. 1  ) 

1760  Rom  =  FNRhoui(T-.i  ) 

■  770  Beta=-2/(Rop+Rom)*(Roc-Rom>/.2 

'780  RETURN  Beta 

•  790  FNEND 

1800  DEF  FNAipha(T) 

1810  Alpha='rNi<u(T)/(FNRhou)(T)-FNCDw(T)  ) 

1820  RETURN  Alpha 

;330  FNEND 

1840  DEF  FNSigma<T> 

1350  Tt=647.3-T 

1360  A-.001*Tt  *2»< . ' 1 60936307/ ( 1  +  .83*Tt>) 

1370  B=.  00 11 21 404638-5. 752805 18E-6*Tt 

'380  C-1.28627465E-8»Tt  2-1 .  14971929E-1  Wt~3 

1890  Siana=A+6+C 

1300  RETURN  Signa  , 

1910  FNEND 
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*•■**    Geometric    'Variables    **~ 

Tur>e   diameter    (D)    =    4.7S3E-03    (m) 

T  _-  -      '  -  -  -  i  '       ;  •    ;  -  r>  1      /  _  > 

I-  /  V  =  \1  CO 

Rotor  radius  <R>   =  0.402  <m) 

Computed  values  are: 

***  .berating  Variables  *** 

Input  variables  are: 

Heat  load             =  50.00  (W> 

Coolant  inlet  temp     =  45.00  (C> 

Coolant  outlet  temp     =  45.72  <C> 

Coolant  mass  flow  rate  =  1.6329E-Q2  (kg/s> 

Fluid  properties  evaluated  at  45.36  (C>  are: 

Specific  neat  <Cp>  =  4224.7  (J/kg.K) 

Viscosity  <Mu>  =  5.866 1E-04  (N.s/m'2) 

Thermal  cond  (if)  =  0.6330  (W/m.K) 

Prandti  number  (Pr>  =  3.884 

Density  <Rho)  =   390.0  (kg/m"3> 

Kinematic  vis  (Nuv)  =  5.9255E-07  (m'2/s> 

Coef  ther  exp  <Beta)  =  4.1266E-04  (1/K) 

***  Calculations  »** 

Preliminary  calculations: 

Reynolds  nunoer  (Re)  =  7.^41E+03 

Friction  factor  (stat)  =  3.402E-02 

Heat  fin*  (Qdp)  =  4.060E+03  (W/m"2) 

M^an  Fluid  vel  <Vm)  =  9.26E-01  (m/s) 

RP«  =  3600.00 

Results: 

Nudb  Nunak        Munt           Nuum            Mus          Nuwm2 

49.52  55.30      111.80        31 .56        44.84      102.34 

Ra  Ko             Ra*Re*pr             M-f                    Nak  1 
7.39E+03     5.16E-01      2.14E+08      1.63E-02        44.90 

■''aRe  ^a                  Ty 

5.50E+07  7.506E-05     0.C26 
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APPENDIX    C 
THERCOSTFHON    ASA1YSIS    (COHPOTSR    PROGRAM    AND    RESULTS) 

The  following  program  was  written  to  give  the  results  of 
the  correlation  for  the  heat  transfer  for  the  rotating, 
ciosed-lcop  thermosyphon  discussed  in  Chapter  III.  It  was 
written  in  HP-Basic  for  the  HP-9826  computer.  Sub-programs 
for  the  graphics  and  the  thermophysical  properties  are  not 
shown.  A  sample  run  is  included  immediately  following  the 
listir.g. 

-Pages      114-117      Program   Listing 

-Paoe        118  SamDle    Calculation 


114 


•  300! 

_  ■ 
1330 
1040 
1050 
iQSO 
1070 
'080 
1090 

*  '00 

mo 

'  '20 
1130 

I  '40 

1150 

-so 

1170 
'180 

;  -90 

1200 
1210 
1220 
:230 
1240 
1250 
1250 
:  270 
1280 
T90 
1300 
1210 
'320 
'330 
'340 
•250 
1360 
'370 
'380 
1390 
'400 
1410 
M20 
1430 
1440 
1450 
*460 
1470 
'480 
1490 
1500 
1510 
'520 
"530 
1540 
1550 
'560 
'570 
'580 


FILE  NAME : 

COM  ZCmin/ 

PRINTER  IS 

BEE° 

PRINT  USING 

D=4.763E-3 

^=.3230 

R=.4r.21 

Kff-0 

i  t  =  2*L  +  ^*;-' 

D?INT  USING 
USING 
USING 
USING 


THSY 

-  e  d  u  r  a  r  v 


PRINT 
PRINT 

PRINT 


Xmin 

1 

2X 

» 


"4X 
"4X 
"4X 
"4X 


Ymin.Sfx,Sfy 

""Default    values 
Tube    diameter 
~ube    length 
Radius    o£_rotor 
r  net  ion   ?  actor 


(m> 
<m) 

(rri) 
for 


benas 


""Tube  diameter 
""Tuoe  lengtn 
""Ckt  Length  <Lt>  = 
""Racius  of  rotor  = 


INPUT 
BEEP 
INPUT 
TF 


.IKE  A 
Ihc-1  THEN 
IF  Id-1  THEN 

BEED 

"ENTER 


OK  TO  ACCEPT  DEFAULT  VALUES 

Q=N>?" 


HARD  COPY  ( l=Y 
PRINTER  IS  701 
280 


1=Y, 

Ihc 


Z.4DE. 
.DD.D." 
,DD.D.' 
,Z.3D,' 

0-N)' 


(m> 

m) 

m) 

m) 


";D 
:L 
;Lt 
:R 


Id 


INPUT 

BEEP 

INPUT 

3EEP 

INPUT 

PRINT 

PPINT 

PRINT 

DRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINTER 

BEEP 

=*RINT 

PRINT 

PRINT 

INDUT 


UBE  DIAMETER  (m>",D 
"ENTER  TUBE  LENGTH  (m)".L 
OF 


"ENTER 

USING 

USING 

USING 

USING 

USING 

■JSING 


RADIUS 
'10X. 

'14-X  , 

'14.X, 
"14X, 
'14X 
■14X, 


ROTOR  (m)".R 

Geometric  Variables  *** 


Tune  diameter  (D>  = 
Tube  lengtn  (L>    = 
C'<  t  lengtr.  (Lt) 
I  /n 

Rotor  raaius  ( R  > 


Z.3DE.' 
.DD.DD.' 

DD.DD.' 
,5D.D":I 
,Z.3D."' 


<m 
(m 

(m 
/D 
<m) 


USING  " 1  OX , ""Computed  values  are 


USING  "2X 

USING  "4X.""1 

USING  "4X.""2 

0p 
IF  0d>1  THEN 
BEEP 
J;  =  0 

PRINTER  IS  1 
INPUT  "LIKE  TO  PLOT 
BEEP 

INPUT  "Nudb-1 .Nunak: 
IF  Okplot-1  THEN 
CALL  Plot(Jj.X.Y.Dt.Tc) 
INPUT  "LIKE  A  LABEL7(1=Y.0=N)".I1 
IF  11=1  THEN  CALL  Label 
END  IF 
END  IF 
Nstep=1 1 
IF  0P=2  THEN 
BEEP 
INPUT  "ENTER  RPM  RANGE  (MIM .MAX)" .Romi .Romn 


SELECT  OPTION: 

Single  point 

Nu  versus  RPM 


(1=Y.0  =  N>r\0kplot 

=2 . Nuwm=3 . Nus=4 . Nul =5" . Ian 
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530 

6i0 
620 

630 

640 

650 

560 

670 

630 

590 

700 

710 

720 

730 

74G 

750 

760 

770 

780 

"30 

300 

870 

320 

330 

840 

350 

860 

370 

'880 

1890 

!900 

910 

920 

330 

840 

350 

960 

970 

980 

990 

2000 

20 10 

2020 

2030 

2040 

2050 

2060 

2070 

2080 

2090 

2100 

2110 

2120 

2130 

2140 

2'50 

2160 

2' 70 

2180 

2'90 


EEE^ 

IWPln  "ENTErt  rfPM  fih  unmp"  Do~ 

END  IF 

BEEP 

IF  Ihc-1  THEN  PRINTER  IS  701 

IF  Inc=0  THEN  PRINTER  IS  1 

DRINT 

USING 


***  Operating  Variables  *** 


input  variables  are 


HEAT  LOAD  (W>".Q 
I4X.""Heat  ioad 

COLD ^SIDE  TEMP  <C>'V 
14X,""Coid  side  teno 


PRINT  USING  "1  OX 
PRINT 

PRINT  USING  "10X 

BEEP 

INPUT  "ENTER 

PRINT  USING  " 

BEE3 

INPUT  "ENTER 

PRINT  USING  ' 

Dt-5 

Th=Tc+Dt 

Ta=^Tc+Tn>* .5 

Dtt-Dt/L 

DRINT 

Cp-FNCouCTa) 

Mu=crNMuu)<  Ta) 

K»FNKu< Ta) 

Pr-FNPrw(Ta) 

Rho=PNRhou<Ta> 

Nuv-Mu/Rho 

Beta=FNBeta<Ta> 

0mega=f\pm*2*PI/60 

Th=Tc+Dt 

Ta-Tc+Dt/2 

PRINTER  IS  1 

PRINT  "TA  -  ".Ta 

Rhoh=FNRnou(Th) 

PRINT  "RhOH=".Rnon 

Rhoc=~NRhou<Tc> 

PRINT  "RHOC=".Rhoc 

Rhom=PNRhou(Ta) 

Drho=Rhoc-Rhon 

IF  Drho<0  THEN  GOTO  3140 

PRINT  "dro=".Drho 

Mu=FNMuw<Ta> 

Cd=FNCpw< Ta) 

Beta=FNBeta(Ta> 

Pr=FNPru<Ta) 

Nuv  =  Mu./Rhom 

PRINT    "MU.CP      =".Mu.Cp 

Mt=0/<CD*Dt) 

Re=4/PI*Mf/(Mu-D) 

Gr=R*<0mega"2>«Beta*<  <D/2 > ' 4 )*Dt t/Nuv"2 

Ra=Gr*Pr 

Prod=Ra»Re»Pr 

Nuwm=.262*Prod~  .  173*48/ 11 

Dtuip  =  G/(K>Nuwm»pI*L> 

PRINT    "DTHP    =".Dtu»p 

Dtuj  =  Dtwp 

Grr  =  R*0mega'2*&eta*D   3»Dtu/Nuv  '2 

Na-Pr/Grr»Lt/R»D/L*Rhom/Rhoh«(Kff + (256/Re*Lt/D > >*Re 

Dtwc-9/<PI*L«K«Nu) 


"".4D.DD.""      <W> :Q 

"",4D.DD.""    (C) ;Tc 
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220  n 

222n 

::3G 

2240 

2250 

2260 

2270 

2280 

2290 

2300 

23:0 

2320 

2230 

2340 

2250 

2360 

2370 

2280 

2390 

2400 

2410 

2420 

2430 

2^40! 

2450! 

2460 

2470! 

24S0! 

2490 

2500! 

25*0! 

2520 

2530' 

2540! 

2550 

2560 

2570 

2580 

2590 

2600 

2610! 

2620! 

2630 

2640 

2650! 

2660! 

2570 

2680 

2690 

2700 

2710 

2?20 

2730 

2740 

2750 

Ta 

2760 

2770 


IF  RBS(Dtiuc-DUi>>.'!G  ~ur'.' 
IF  Rpra>3000  THEN 
n  -  =  n +  *  q 

OU  1  u    t  _.  \J  w 

END  If 

IF  Rpm>!000  THEN 

Dt-Dt+1 

GOTO  '900 

END  IF 

IF  Rpm>500  THEN 

Dt-Dt+1 .5 

GOTO  1900 

END  IF 

IF  RDm>=300  THEN 

Dt=Dt+4 

GOTO  !900 

END  IF 

END  IF 

V««4«Mf/<Rhom*PI»D"2) 

Ro=Vm/(0mega*D> 

Jay=0mega*D "2/Nuv   !From  STEPH-Rotat lonal  Re 

Dtt-<Th-Tc)/L 

Gr»R*< Omega ~2>*6eta*<  (D/2)  '4  >«Dtt/Nuv  '2 

Ra=Gr»Pr 

Laminar    Value 

Nulam=48/1  1 

Turbulent  Value  from  Di t tus-Boe 1 ter 
Nudb= . 023*Re  .  8*Pr'  .4 

Stepnenson  Correlation 

Nus-<Pr    .4/ .367>-.007l*Re    .88*Jay~ .023    !  Air  ,Tu.rb-EQN2.23-corrected   for   uial 

Naicavama  Correlation    *or    ianinar    conditions 

J2-J*2 

Sq-2/1  !•<  :-rSQR<  1+77/4»(1/Pr"2))l 

Cf  =  1  -  .486-( (3»Sq-1 )'  .4/<Sq«<Sq*Pr*S0R(5)+2)>)*J2/(Ra»Re)    .6 

Nunax 1=48/! 1». 191 /Sq*<3*Sq-l  >  ' . 2*<Ra-Re>  '  .2*1 /( i+< 1 /10~Sq*Pr > ) 

1us=ttu    ! Temp   Vaiue   MusiMu 

Prod=Ra*Re*Dr 

Woods-lorr is  Correlation  for  Laminar 
Nuwm=.262»Prod' .173-48/1 1  !  Eqn  2.22-WM3 
Jl=Jay/8 

Woods-Morris  Correlation  for  Radial  Pipes 

Nuwm2=.015*Re' . 78*J1 ' .25  !from  um2 

IF  lhc=0  THEN 

PRINTER  IS  1 

ELSE 

PRINTER  IS  701 

END  IF 

IF  Jj=0  THEN 

IF  Jj  =  0  THEN  PRINT  USING  "1  OX . ""Resul ts :  

PRINT  USING  "10X.""Fiuid  properties  evaluated  at 


DDD.DD,""  (C)  are 


PRINT  USING  "14X.  ""Specif ic  heat  <Cp> 
PRINT  USING  "14X. ""Viscosity  <Mu) 


"",4D.DE.""<J/kg.K> :Cp 

■•■'  .Z.4DE.""  (N.s/nT  2) :Mu 
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278( 
2790 
2800 


2830 
2340 
2S50 
2SG0 
2870 
2880 
28^0 
290C 
2910 
2920 
2930 
2^0 
2950 
2S60 
2970 
2980 
2890 
3000 
3010 
3020 
3030 
3040 
3050 
30G0 
3070 
3030 
3090 
3100 
3110 
3*2C 
3130 
3140 
3150 
3.!  50 
3170 
2180 
3190 


PRINT 
PRINT 
PRINT 


USING 
USING 

USING 


USInG 

USING 
USING 
USING 
USING 
USING 


"14X 
"14X 

"14X 

"I  ■ 

'  I4X 


" "P rand  1 1 

""[■ens  l  ty 
■\  x  nema  t  .  ^ 

Lee*  Therm 


cord  ( k > 
number  (! 

- 


_..pi  oera> 
"1 4X  ,  ""Reynolds  number  (Re> 
4X.'",s1ean  fluid  vel  (Vm) 
'"Mass  fiou  rate 
""Temperature  (hotside) 
""RPH 


"MX 
"14X 
"14X 


PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

END  IF 

PRINT 

"    3lcpiot=0   OR    (Olcpiot-1      ME 

PRINT  USING  "10X.""    Nudb 

=>RINT  USING  "13X.S<3D.DD.2X>' 

PRINT  USING  "14X.""   Ra 

PRINT  USING  "14X.4(Z.2DE.2X> 

END  IF 

IF  Qkplot-1  THEN 

Ian-1  THEN  Y=Nudb 

Ian  =  2  THEN  Y  =  Nunaicl 

Ian=3  THEN  Y=Nuun 

Ian=4  THEN  Y=Nus 

Ian  =  5  THEN  Y=( 48/11 > 
X=Rpm 
J  t  =  J  j  + 1 

CALL  Plot(Jj.X.Y.Dt.Tc) 
DRINTER  IS  1 
END  IF 

IF  0p=2  THEN 
Dt-5 

Rom  =  RDn*1 0'  . 1  - 
IF  Rpm>Rpmh  THEN  3130 
GOTO  TS9c 

END  IF 

BEE^ 

INPUT  "ANOTHER  RUN  ( 1 -Y  ,0=N>?" . Ir 

PRINT  "PIT 

IF  Ir-1  THEN  1360 

INPUT  "WANT  TO  LftBEL?<1-Y.0-N)".I 

IF  11  =  '  "THEN  CALL  ^aDel 

END 


Z.4D."" 

Z.3DE":P 

4D.D."" 

_     .    —1  L,   _     t 

L  ■  MUC  .  " 

Z.3DE":R 

Z.2DE."" 

Z.2DE."" 

DDD.DD." 

DDDDD.  ." 

< W/m.k » :< 

T 

r*  3/ m  3  I """ : Rho 

(  1  /K)       ' : Beta 

e 

(m/s) :Vm 

(kg/s) :Mf 

"(C) :Th 

:  Rom 


lstep<1  !  >  THEN 

Nun  air  1    Nuwm      Nus 
:  Nudb  .  Nunak  1  ,  Nuum  .  Nus 

Ro     Ra*Re»Pr     Mf 
■.D.DD":Ra.Ro.Prod.Hf 


IF 
IF 
IP 

IF 

IF 
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*■*■"    is  :  '<  e  t  r  1  c    >  a :  .  _l_  ,.  c  -■  « *  * 

ToDe   diameter    lU)    =  3./b3E-03    (m) 
Tube    length    (L)         =         . 32    <m) 
Oct    length    (Lt>         =      2.45    <n> 
L/D                                 =        172.8 

Rotor    radius    (R)      =  0.402    (m> 

Ccnouted   value?    are: 

***    Operating    !-'ar*ar>ies    *** 

Input    variables    are: 

Heat    load  =        50.00      <w> 

Cold    siae    temp  =        45.00    (C) 

Resul ts : 

fluid   prooerties    evaluated  at      £2.75    (C>    are: 

Specific    neat    (Cp)  =    4? i 0 . 1 E+00 ( J/kg . K > 

Viscosity    <Mu>  =    4.4172E-04    (N.s/m   '. 

Thermal    cond    (k)  =    0.8404    (W/m.k) 

Prandtl    numDer    <Pr>  =    2.302E+00 

Density    (Rno>  =      383.1    9kg/nT3> 

Kinematic    vis    (Nuv>  =    4.4991E-07    (M"2/S> 

Coef    therm    exp(Beta)  -    5.3940E-04    (1/K) 

Reynolds    nomoer    (Re)  =    I.978E+02 

Mean    fluid   vel    <Vm>  =    1 .87E-02    (m/s) 

lass    flow    rate  =    3.27E-0*    (kg/s) 

Temperature    (hotside)  =      80.50(C) 

RPM  =      3600. 

Nudb  Nun  at  1         Nuum  Nus 

2 .42        5 i .30        34 .4!  I .65 

Ra  Ro  Ra»Re*&r  Mf 

S.13E+05      1.04E-02      3.52E+08      3.27E-04 
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APPENDIX   J> 
CLOSED    TSO-FEASE   THE E20SYPH ON    ANALYSIS     (COMPUTER    PROGRAM    AND 

RESULTS) 

The  following  program  was  written  to  give  the  results  of 
the  correlations  for  the  sonic,  boiling,  and  condensing 
limits      as    discussed      in     Chapter  IV.  It      was    written     in 

HP-Basic  for  the  HP-9826  computer.  Sub-programs  for  the 
graphics  and  the  t hermophysical  properties  are  not  shown.  A 
sample    run    is  included   immediately   following   the   listing. 

-Fages      120-12U      Frogram   Listing 

-Pages      125-127      Sample    Calculation 
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'   ■ 

'  )4  ■ 
r  050 

■  -  ] 

1090 

"  '00 

'  •  -  ' 

•  2C 

•  '  30 
1140 

•  !50 
1160 

I  !70 

II  SO 

i  :90 

1200 

• :  • : 

1 230 
■240 
1250 
1260 

•270 

• :  ■ 
;22o" 

1300 

1320 
K330 
1340 

;  ?50 

1360 
1370 

H3  0 
1390 
1400 

14J0 

1420 
1430 

14*0 

1450 

1^60 
1470 
I  480 
1490 
1500 
1510 
1520 
'530 
1540 
1550 
1560 

•  =  70 
'530 


.  n  I  - "  1 I  ;  M.  a  r  c  h  5  .  '  - 
-AS  '  RE  '  "-TI  ■  '"-.-  -  -  *  ;  ~  "•■ : 
2170.CALCS:  4000. PLQ"^  5250. FNS:  61 00 .LABEL: 

riu   Cmn/  <nin.  Ynin  .Sf  x  .of  y 

PRINTER  IS  1 

3EE3 

DRINT  USING  '^X.'^'Defauit  values:  

D=4.763E-2     !  Tone  diameter    <m> 

!_?=. 3230/2    !  Tube  lengtn(EvapMn) 

_  c  = . 5  *L e 

_  t  -  L  e  ~"_  c 

R =  •  - !"  2  '        !  R  ac:  i  us  o  f  r  o  •:  or  <  ,t>  > 

PRINT  JSING  "AX.""Ti±De  diameter    =  "".Z.4DE.""  Cm)"" 

PRIN1  JSIMG  "4X.""Tube  Iength<Evap)=  "".DD.D.""  <m> 

PRINT  USING  "4X.    ube  Lenath<Cond)=  "".DD.D.""  <m> 

PRINT  USING  "4X.""3adius  or'    rotor   =  "".Z.ID.  ""  <-i> 

BEE3 

INPUT  "G:<  70  ACCEPT  DEFAULT  VALUES  M  =Y  .0=N>?" .  Id 

3EE3 

rMPUT  »!_XKE  J  HARD  :2DV  (•  =  y  .'  =  '-!)  °".  !hc 

Ir  Ihc=1  THEN 

-"INTER  is  70! 

ELSE 

PRINTER  IS  I 

END  :r 

[F  id=l  "HEN  1360 

INPUT  "ENTER  "USE  DIAMETER  <m>".D 

2EE- 

:  PUT  "ENTER  TUBE  LENGTH<E>/ap>  (m)".Le 

3EE= 

INP'JT    "EvlrER    T'JBE    _ENGTH<Ccna>    <n>".Lc 

3EE° 

INPUT  "'ENTER  RADIUS  OF  ROTOR  Cm>".R 

PRINT  JSING  ";0X.""»»*  Geometric  Variables  *»*  

PRINT -USING  "14X.""Tube  diameter  (D)  =  "".Z.3DE.""  Cm) 

PRINT  USING  "14X.""Tube  ienath<Evap>  (Le^   =  "".Z.3D." 

Dt?lNT  USING  "14X.""Tube  lenatn<Cond>  <  Lc )   =  "".Z.3D." 

PRIN~  'JSING  "14X. ""Rotor  radius  (R>        =  "".Z.3D." 

PRINT 

PRINTER  IS  1 

3EE3 

PRINT  JSIMG  "2X .  "'"SELECT  OPTION: 

PRINT  USING  "4X.""'   u  versus  %Fill 

oc?INT  USING  "4X.""2   77*??????'???? 

INPUT  Dp 
IF  0d>!  THEN 
3EE° 
PRINTER  IS  I 

PRINT  USING  "4X.""N0T  IMPLEMENTED  YET 

GOtq  1  iiAij 
ELSE 
BEE- 
G-0 

Jj-0 

P  -  r  ■  • 

l"lP'jT  "LIKE  TO  PLOT  (l=Y.0=N)?".0kDiot 


":D 

(m)' 

(m)' 

<m)' 
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-  C  _  ■ 

Z-     Dk  o  i  o  t  =   "HEN 

' 

iLL  :  Lc  :  ■  .;  .Jc.Jd.  <.  f  .Type.  :x) 

R^ra 

INF J7  "1-i JLIl  . 2=DRSn-0ASH . 3=DQT-DC  f"  ,  Type 

-:- : 

END  "r 

IS50 

END  IF 

1660 

\- teo=100 

1570 

BEEP 

*660 

INPUT  '"ENTER  OPTION  ( 1=S.2=E.3=3.4=C>".  Io 

1690 

Ic  "c  =  2  THE'. 

;  700 

3EED 

!7TG 

INPUT  ,",=WALLlS.2=JASTER".Ie 

_ 

END  IF 

73C 

Zr  :-.:=•  "HEM  PRINTER  13  701 

1740 

IF  lhc=0  THEN  PRINTER  13  1 

,T50 

BEEP 

•73o 

IMF  JT  "ENTER  "sat  CO". Tsat 

•7^0 

PRINT  USING  "14X.""Tsat<C> 

"",DDD.D":Tsat 

•  730 

BEEP 

1790 

INPUT  "ENTER  RP^'-Rom 

'300 

PRINT  JSING  ■••ax.,,,,R='M 

"" . 4D" :Rcm 

1  8  1  0 

3EE- 

'320 

INPUT  "ENTER  ;0ND.  WALL  TEMPERATURE" . Tu 

1830 

PRINT  JSING  "14X.""T  ua!  . 

"".DDD.D":Tuj 

:3d0 

At-PI»D'2/4 

1850 

Tsatk=Tsat+273.l5 

'360 

Cp-h NCou<  sat) 

•37^ 

Mul=FNMotij<  (Tsat  +  Tu)  72  ! 

!?30 

P1ij.v=p'NMuv  <  Tsat ) 

1890 

K=FNKu< Tsat  > 

1900 

Pr-FNPrw< T;at) 

1910 

Rhoi=FNRhow<Tsat> 

1920 

Hfg-P"NHfg<Tsat) 

1930 

Beta-FNBeta< i sat > 

'OUQ 

Sigma=1~NSigma(  sat) 

1950 

Rhov=" /FNVvst < Tsat ) 

i  350 

Mu.v  =  Ma.//Rhov 

1970 

Rsh-1.327 

1980 

Cv-Cp/Rsh 

1990 

Rgc  =  «61  -52   !  CJ/kg.K) 

2000 

Clmega~Rpm»2*PI/6G 

2010 

Accel=Omega  2*3 

2020 

IF  Jj=0  THEN 

2030 

PRINT 

2040 

PRINT  USING  "1  OX , ""F laid  properties  evaluate* 

3  at  ""  .DDD.DD.  "" 

(C)  are: 

2050 

PRINT  USING  "I4X. ""Specif ic  heat  <Cp)     -  ' 

'••.■^D.D.,•'•  (J/kq.K) 

;Cp 

2060 

=>d:,i"  'JSING  "14X  ,"MSa-io  of  Specific  Heats-  ' 

"\D.DDD":Rsh 

2070 

PRINT  USING  "14X. Viscosity  (Ma;         -  ' 

■".2.4DE.""  (N.s/m 

2030 

PRINT  USING  "14X. ""Thermal  cond  <k> 

'".Z.4D.""  (H/m.K)' 

2030 

PRINT  USING  "14X.""Prandti  numoer  (Pr)    =  ' 

"\Z.3D":Pr 

2100 

PRINT  USING  "14X. ""Enthalpy  (Hfg)         =  ' 

•".Z.4DE":Hfa 

2"0 

PRINT  USING  "i4X.""Density  (Rhoi)         =  ' 

"\4D.D.""  (kg/ra  3) 

;Rhoi 

2'20 

PRINT  USING  "14X. ""Density  (Rhov) 

"\Z.4D.""  (ka/m  3' 

"""  :Rhov 

2*30 

PRINT  USING  "14X.""Coet  tner  exp  'Beta)    =  ' 

'".Z.4DE.""  (1/K)4" 

"•.Beta 

2'40 

PRINT  USING  "•4X.,,,,0meaa 

'".4D.D.""  (1/s) 

' :  Omega 

2'50 

PRINT  USING  "!4X. ""Surf ace  Tension 

■".Z.4DE.""   (N/m)' 

"" : S  igma 

2160 

END  :- 

2170! 

?  1 30 ! 

pre   =  'cj 
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2"  Q   • 

22C  ;   i=Rh  :  .-_:-->-  -_r     ■   ;;  *Q  I'  '  J 

2210   flv=(M/(Lt*«  Rhov-Rho  [)>>-<  (Rhoi  /"(Rhov-Rhoi  >  >*At> 

2220! 

2230!  SONIC  OAlCLL- :  ION 

2240   IF  Io=1  THEN 

2250   IF  J •=:  THEM 

2250   PRINT 

2270  PRINT  JSIMG  "tOX. ""Sonic  Limit  Calculations: 

2230   PRIM"  USING  "10X.""*  Fill    Tctai  M€at 

2290   END  IF 

Om=R  hov  «*ifg*<  1/ <  Lt*(  Rhov-Rho  1)  >-Rho  I  /(Rhov-Rho  1  )*At  >*SGR<  Rsh*Rgc*T=a: 
2310   Mf=Qm/Hfg 
232   X=?cf 
2230   r=0m 

2340   Ir  3n<=0  THEN  GOTO  2^00 
2350   IF  0kpiot=1  THEN 
23S0   CALL  Dio:<o;.~c.-c.X.v.Tyoe,Cx> 
2370  ELSE 

2330  PRINT  JSIMG  "!2X. 3D.8X .Z.3DE" :Pcf .Qm 
2390   END  IF 

_   ]   IF  Pcf>=!00  THEN  GOTO  3730 
241C   Ir  0kplot=1  THEN 
2420  Pcf=Pc^*l0".02 
2430  ELSE 
2440   Pcf=Pcf+l 
2450  Jj-jj+1 
2460  END  IF 
2470   SOTC  220C 
2*90   END  IF 
2490! 

2500!  ENTRAINMENT  :ALCULATIGNS 
25: 0   IF  Io-2  THEN 
2520! 

2530   IF  Jj-0  "HEM 
25Z.0   PRINT 

2530   PRINT  USING  "!0X. ""Entraining  Limit:  Ciaculations: 

2550   PRINT  USING  "10X.""JJ  riii     Total  Heat 

2570   END  IF 

2530'    wALl.IS   CORRELATION    J*=.25 

2590      IF    Ie-1    THE'. 

2S00   IF  Av«0  THEM  GOTO  2090 

2610  Jv=.25/Rhov  .5»(0mega'2*R»D»(Rhol-Rhov>)".5 

2F20   Qm  =  Rhov-»  Jv*At*Hf  g 

2530  END  IF 

2640! 

2550!  JASTER  CORRELATION 

2650   IF  Ie=2  THEN 

2570   IF  Pcf-0  "HEM  Pcf-I 

2530  T1*(1/(1-(Pcf/100>))-1 

2530!  T1-<1/(Pcf/100>>-1 

2700   T2=x1*(Rhol/Rhov)-<2/3> 

27to   IF  "!<0  THEM  2900 

2720  Xx2-<Mil1/Muv>  *.25*T1  I  .75*<Rhol/Rhov>  s<  <2/3)*1  .75-1 ) 

2730  Xx-1/(1+T2) 

2740*  <x=.5 

2750!    <x2=<r4ul/Muv)  "  .25~Rhov/Rhol 

2750      IF    G=0    THEN   G=1 

2770      Rel=<?-Xx)*<G*D/Mul> 

2730      Rev=Xx*G*D/Muv 

2790      Fv=.079/Rev    .25 
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2S    :  IF    Re<=1250    THEN   Dd='SGR(Rei/2> 

28!Q  IF    3e>"25'G      HEN   Dd=.G504*<Rei>     .375 

2820  Num=1Q*SQR<2>*<Rhoi*Rhov>  ! .5*Mui*Dd*0mega  2~- 

)  Den=(  !  ^<<2  '  .25)  S*RhoI*Xx  "--•<  (3  2> 

2S^u  Ge=<Num/Den>  (1,  5) 

2S5G  IF  ABS<<G-Gc>/Gc>>.0!  THEN 

2860  G=<G+Gc>/2 

2370  GOTO  ?7£Q 

2830  END  IF"' 

2880  Qm=Xx*G*At*Hfg 

2900  END  Ir 

2910  X-Pcf 

2820  Y=0n 

2830  IF  3n<  =  0  THEN  GOTO  2890 

2940  IF  Okpiot=1  THEN 

2950  CALL  PIot< J j , Jc. Jd.X. Y.Type .Cx) 

2'?i:  ELSE 

2970  "-INT  USING  "12X.3D.3X,  Z.4DE.4X.Z.3DE"  :?cf  .Qm.G 

29P-C  END  "r 

2990  IF  Pcf>=f00  THEM  GOTO  2730 

3000  IF  GkDiot=1  THEN 

3010  Pcf=Pcf*IO\02 

3020  ELSE 

3030  Pcf=Pcf+l0 

30^0  END  IF 

3050  .•=.--• 

;:-:o  :r  ie=2  them  ;g~2  :c,o0 

3070  IF  Ie=1  THEN  GOTO  2200 

3080  END  l- 

3090! 

2 '00!  BOILING  CALCULATIONS 

3110  IF  Io=3  THEN 

2-20  IF  Jj=0  'HEN 

3130  PRINT 

3J«0  PRINT  USING  ,,10X.""3oi  ling  Lim'it  Calculations: 

3150  D»INT  USING  ■"  OX.  "'"/.    -  l  I  L     Totai  heat       QDP 

2 '50  END  IF 

2 ' 70  IF  Pcf>!00  THEM  GOTO  3720 

3180  Ab  =  D*'_e*ACS( 1  -Pcf  ^5Q) 

3190  Qm=Ab*.  ;2-^K.ov  .5*Hfg*<0mega*2*R*<Rhol-Rhov)*Sigma)  .25 

3200  X-Pcf 

3210  r-Qn 

2220  IF  Qra<-0  THEM  2230 

2230  IF  Olcplot»1  THEN 

3240  CALL  PloUJj.Jc.Jd.X.Y.TyDe.Cx) 

3250  ELSE 

32S0  PRINT  JSING  "12X.3D.SX.Z .4DE.3X,Z.4DE":Pcf .Qm.Qm/Ab 

3270  END  IF 

2230  IF  PcfMOO  T^EN  3-OTO  3730 

3290  Jj=Jj+1 

2200  IF  2lcDiot='  THEN 

3310   Pcf=Pcf*l0~.02 

3320   ELSE 

2330   Pcf=Pcf+1 

3340  END  IF 

3350  GOTO  2200 

3360  END  IF 

3370' 

1330'  CONDENSATION  CALCULATIONS 

2290   IF  Io-4  *HE'i 
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.-•  "i  IF  !j=(  "HEN 

;_•  PRINT 

3420  PRINT  JSING  '  *  OX  .""Cortdensina  Limit  Caiculat  ions  : 

3430  PRINT  USING  "'OX. ""7,  "ill    "Total  Heat 

3440  END  IE 

3450  Dt=Tsat-Tu 

3^60  HfgD=rifg*( 1 »■  .68*<Cp*Dt/Hfg) > 

3470  Tag=(Rhol*(Rhoi-Rhov)*Qmega  2*R*Hf  ap*k'"3/(D*Mi±i*Dt ) )  .25 

:  - : ;  deg 

3490  Phi=!30-ACS< ! -<Pcf/50>) 

3500  IF  Phi>=0  AND  Phi <=1 1 0  THEN  F=FNF1(Phi) 

3510  Ir  Phi> 1 10  AND  Phi<=!50  THEM  F-FNF2<Phi) 

3520  IF  Dhi>i50  AND  Phi<=;30  THEN  F=FNF3<Phi) 

3530  IF  °hi>!S0  THEN  PRINT  "ERROR .PHI> I  30 ! ! ! " 

35-0  RAD 

3550  Qm=Av*Dt*F*Tag 

3560  X=Pcf 

3570  Y-Qn 

3580  Ir  3m<=0  THEN  GOTO  3700 

3590  :~    3kpiot=1  THEN 

3600  CALL  PIot( Jj . Jc.  Jd.X. Y.Type.Cx) 

3610  .  E_3E 

3620  PRINT  USING  ,,!2X.3D.8X.3D.DM:Pcf  .Qm 

3630  END  Ir 

3S40  IF  0kpiot=!  'HEN 

3650  Pcf=Pcf*10-.02 

3660  ELSE 

3670  Pcf-Pcf+1 

3630  END  Ic 

36 30  Jj=Jj+1 

3700  IF  Pcf>=!00  THEM  GOTO  37  30 

3710  GOTO  2200 

3720  END  IF 

3730  3-rD 

3740  INPUT  "ANOTHER  RUN  U=Y.0=N>?",Ir 

3750  IF  Ir-1  THEN  120G 

3750  3EEP 

37^0  INPUT  "WANT  _;3  LABELS !-Y.0-N>*MI 

3730   IF  11=1  THEN  CALL  Lacei 

3730  END 
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Jeoraetnc  /ariaoies  *~* 
Tube  diameter  (D>  =  4.763E 

•_re  length (E^ap.)  <!_-9> 
Tube  Length(Cond)  <Lc)   =    .1 
SoTor  radius  (R)   =  0.402  (m) 


'A\     ■  n) 

21     <  m  > 


'ia->  : 


80 
2600 


:    properties    e  v/  a  i  ua  t  ea 
Specific    heat    (Cp) 


-'i 


=  44c 


-•  J  '  - )  are 
3  (  J  /k  q  .  K .) 


Ratio  of  Specific  Heat< 

Viscosity  (Mu> 

The rmai  cone  < k > 

Prandti  numDer  (Pr) 

Enthaioy  (Hfg) 

Density  <Rnol> 

Densi by  (Rhov> 

Coef  ther  exp  (Beta> 

Omega 

Surface  i  ens  ion 


i  ,327 

4.4336E-04  (N.s/m 

0.6638  <W/m.K> 

2.2'- 
2.3057E+06 

97!  .3    (kg/m   3> 
.2    (kg/m   3) 
6.4578E-04    (1/K> 

377.0    H/S) 
7.7043E-02      (N/m) 


2) 
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Lonaensina 

Limit  Calculation 

'/.    "ill 

Total  Heat 

1 

56.  1 

2 

64.7 

3 

h3.4 

4 

82.2 

0 

61.1 

6 

60  .  " 

7 

58 . 1 

8 

53.8 

3 

55.5 

10 

54.4 

4  « 

CO   ") 

i  1 

.;  O  .  <L 

12 

52.1 

13 

51  .0 

14 

49.9 

15 

48.  9 

16 

47.3 

!7 

46.9 

18 

45.9 

1  9 

45.0 

20 

44.0 

21 

43.1 

22 

42.2 

oo 

41  .3 

24 

40.5 

25 

33.6 

26 

38.8 

27 

37.9 

23 

37  . " 

?Q 

36 .  ? 

30 

35.5 

:  • 

34.7 

0£ 

3°. .  9 

0? 

33."  4 

34 

32.5 

35 

31  .7 

36 

30.8 

37 

30.0 

38 

29 . 2 

39 

23.4 

40 

27.5 

41 

26.9 

42 

26.1 

43 

25.4 

44 

24.6 

45 

23.9 

46 

93  ,2 

47 

22  .'5 

48 

21  .8 
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43 

2!  .  ! 

50 

20.^ 

51 

13.3 

52 

19.1 

53 

'c> . 5 

54 

17.9 

55 

17.2 

55 

16.6 

57 

16.0 

58 

15.4 

59 

14.9 

SO 

•4.3 

6! 

13.7 

52 

13.2 

S3 

12.6 

64 

12.1 

65 

11  .6 

56 

11.1 

67 

10.5 

58 

10.1 

S3 

9.6 

70 

9.1 

71 

8.6 

72 

8.2 

7? 

■  J 

7.7 

74 

7.3 

75 

5.9 

76 

6.5 

77 

6.0 

73 

5.6 

79 

5.3 

80 

4.9 

31 

4.5 

82 

4  .2 

33 

:  .  o 

84 

3.5 

35 

~i   ^ 

■!•   .  d 

86 

2.8 

37 

2.5 

38 

o  n 

39 

90 

1  .7 

91 

•  .5 

92 

*  .  2 

o^ 

1  .0 

-/   sJ 

94 

.8 

35 

.6 

96 

.4 

97 

.  3 

98 

9 

99 

.7 
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